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EXPERIMENTAL AND PREDICTED LONGITUDINAL AND LATERAL-DIRECTIONAL RESPONSE 
CHARACTERISTICS OF A LARGE FLEXIBLE 35° SWEPT-WING AIRPLANE 

AT AN ALTITUDE OF 35,000 FEET 1 

By Henry A. Cole, Jr., Stuart C. Brown, and Euclid C. Holleman 


SUMMARY 

Measured and predicted dynamic response characteristics of 
a large flexible swept-wing airplane to control surface inputs 
are presented for flight conditions o f 0.6 to 0.85 Mach number 
at an altitude of 35 ,000 feet. The report is divided into two 
parts. The first part deals with the response of the airplane to 
elevator control inputs with principal responses contained in a 
band of frequencies including the longitudinal short-period 
mode and several symmetrical structural modes. The second 
part deals with the response of the airplane to aileron and 
rudder control inputs with principal responses contained in a 
band of frequencies including the Dutch roll mode , the rolling 
mode , and three antisymmetrical structural modes. 

The measured dynamic response characteristics , excited by 
pulse-type control inputs , are presented in frequency-response 
form for response quantities at the center of gravity , several 
locations on the wing , and the tail. The low-frequency part of 
the response {frequencies below natural frequencies of structural 
modes) is also presented in transfer function form and com- 
parisons are made with transfer functions predicted by equa- 
tions of motion which include first-order effects of flexibility. 
Reasonably good agreement is obtained between experiment and 
theory for frequencies below the natural frequencies of the 
structural modes where the assumptions of the equations are 
valid. Effects of flexibility are shown by comparisons of pre- 
dicted rigid airplane response with the measured and predicted 
flexible airplane response. 

A method for determination of aerodynamic lift and moment 
of a flexible wing through use of aerodynamic and structural 
influence coefficients is presented. Also , relationships between 
transfer function coefficients and stability derivatives are shown, 
which in some cases allow stability derivatives to be approxi- 
mated from flight data. 

INTRODUCTION 

The desire to increase the range and speed of large air- 
planes led to configurations with wings of high aspect ratio, 
thin airfoils, and fuselages of high fineness ratio. All of 
these factors tend to increase the flexibility of the structure, 
and the associated aeroelastic effects become of greater 
importance in problems of static and dynamic stability and 
control. The dynamic effects are especially important when 
the airplane is equipped with automatic control because 
structural modes may introduce system instabilities which 


would not occur in a rigid airplane. The prediction of 
these aeroelastic effects is important for rational design of 
the airplane and its control system. 

In view of the above problems, the NACA conducted 
flight tests on a large flexible 35° swept-wing airplane over 
a wide range of flight conditions. The aims of this program 
were to document and analyze the airplane response to con- 
trol surface motions and, through comparisons between 
measured and predicted response characteristics, to establish 
simple but adequate methods of prediction for flexible air- 
planes. In the present report the longitudinal and lateral- 
directional responses to control surface motions are docu- 
mented at an altitude of 35,000 feet over a Mach number 
range of 0.6 to 0.85. Response quantities at a number of 
locations on the wing and fuselage are shown so that the 
response of the entire airplane is fairly well defined. Also, 
a simple method is developed for the prediction of the 
response for inputs which are primarily influenced by the 
low-frequency modes, the longitudinal short period, the 
Dutch roll, and rolling modes. 

With regard to measuring and analytical techniques, the 
frequency response is a convenient way to express the 
dynamic-response characteristics in a standard manner 
which is independent of the particular input from the control 
surface. Transient response to pilot-applied pulse inputs 
were measured and transformed to frequency response by 
means of the Fourier integral, a technique described in 
references 1, 2, 3, and 4. Although these frequency- 
response data completely define the response characteristics, 
operational expressions relating the output response to an 
input disturbance, known as transfer functions, are a more 
useful form for detailed analysis or for the synthesis of 
automatic control systems. These transfer functions may 
be approximated from the frequency response by a curve- 
fitting procedure such as described in reference 5 or by 
matching time histories on an electronic analog as described 
in reference 3. Comparable predicted transfer functions 
may be derived from the equations of motion including 
rigid body and structural degrees of freedom. 

Simplifications in these transfer functions appeared justi- 
fied when the response of interest wtis limited to the band 
of frequencies bekw the structural modes. Consequently, 
theoretical transfer functions were derived from the rigid 
airplane stability equations with first-order effects of flexi- 


1 Supersedes NACA RM A54H09 by Henry A. Cole, Jr., Stuart C. Brown, and Euclid C. Holleman, 1954, and TN 3874 by Stuart C. Brown and Euclid C. Holleman, 1956. 
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bility included. Coefficients were calculated for both the 
rigid and flexible airplane for comparison with coefficients 
evaluated from flight-test data. Simplified relationships 
between transfer-function coefficients and aerodynamic de- 
rivatives are used in some cases to obtain aerodynamic 
derivatives for comparison with estimated derivatives. 

The data used in this report were obtained from flight 
tests conducted at the NACA High-Speed Flight Station 
and the analysis and reduction of data were a cooperative 
effort of HSFS and Ames Aeronautical Laboratory. 

TEST EQUIPMENT 

The test airplane was a Boeing B-47A with General Elec- 
tric J-47-GE-23 turbojets and with wing vortex generators 
as shown in figure 1. Physical characteristics are listed in 
table I. The airplane was fitted with a nose boom for meas- 
uring airspeed, altitude, angle of attack, and angle of side- 
slip, and an optigraph for measuring the movements of 
target lights on the wing and tail (figs. 1 and 2). Left and 
right aileron, left and right elevator, and rudder deflections 
were each measured by three NACA resistance-type control- 
deflection indicators located at root, midsemispan, and tip. 
The outputs were recorded on Weston 12-channel and Con- 
solidated 18-channel oscillographs. The pitch, roll, and 
yaw rates at center of gravity and roll and yaw rates at the 
tail were measured by self-recording, magnetically damped 
NACA turn meters, the accelerations at center of gravity 
and tail by NACA air-damped accelerometers, and the accel- 
erations on the wing by Statham linear accelerometers. The 
instruments were aligned with respect to the body reference 
line and locations of the pertinent instruments are shown in 
figure 2. 

TEST PROCEDURE 

The flight -test conditions covered Mach numbers from 
0.6 to 0.85 at an altitude of 35,000 feet. For longitudinal 
maneuvers the center of gravity varied from 12 to 30 percent 
of the mean aerodynamic chord and the gross weight from 
110,000 to 120,000 pounds (table II). For lateral-direc- 
tional maneuvers the center of gravity was approximately 
21 percent of the mean aerodynamic chord and the gross 
weight approximately 115,000 pounds. 



TABLE I.— PHYSICAL CHARACTERISTICS OF THE 
TEST AIRPLANE 


Wing 

Area, sq ft 1, 428 

Span, ft 116 

Aspect ratio _ _ 9. 43 

Taper ratio 0. 42 

Mean aerodynamic chord, ft 13. 0 

Sweepback of quarter-chord line, deg 35 

Airfoil thickness ratio (parallel to body center 

line), percent 12 

Dihedral, deg 0 

Ailerons 

Area aft of hinge center line (each), sq ft 53. 8 

Aileron span to wing span ratio_ 0. 405 

Average aileron chord to wing chord ratio 0. 26 

Vertical tail 

Area (including dorsal), sq ft_ __ 227 

Span, ft 18. 7 

Aspect ratio 1. 54 

Taper ratio 0. 34 

Mean aerodynamic chord, ft 13. 02 

Sweepback of quarter-chord line, deg_ 34. 9 

Airfoil thickness ratio (parallel to body center 

line), percent _ __ 10 

Distance )4-chord wing M. A. C. to ^-cliord 

vertical tail M. A. C., ft. 46. 5 

Rudder 

Area aft of hinge center line, sq ft __ . 51. 2 

Average rudder chord to vertical tail chord ratio. 0. 30 

Average gross weight-, lb_ ... 115,000 

Moment of inertia about the X axis, slug-ft 2 __ 1, 074, 000 

Moment of inertia about the Z axis, slug-ft 2 .. . 2, 306, 000 

Inclination of principal longitudinal axis with respect 

to fuselage reference line, deg _ —2. 6 

Horizontal tail 

Area, sq ft 2r>s 

Span, ft 33 

Aspect ratio 4.06 

Taper ratio 0. 123 

Mean aerodynamic chord, ft 8. 57 

Sweepback of quarter-chord line, deg. 32. 9 

Airfoil thickness ratio (parallel to body center 

line), percent 10 

Distance }£-chord wing M. A. C. to ^-chord 

horizontal t ail M . A. ( ft 16. 6 

Elevator 

Area aft of hinge center line, sq ft 68. 5 

Average elevator chord to horizontal tail chord 

rat io 0. 30 


Before each maneuver, the airplane was trimmed at the 
desired speed and altitude. Then the pilot applied a cpiick 
elevator, aileron, or rudder pulse and held controls fixed 
until the transient motion damped out or, in the case of 
lightly damped lateral oscillations (yaw damper off), the 
airplane was allowed to oscillate for about 25 seconds. 
Typical time histories are shown in figures 3, 4, and 5. 

The longitudinal moment of inertia was measured in 
ground oscillation tests in which the airplane was supported 
on knife edges and a spring (ref. 6). Other inertia charac- 
teristics were estimated. 

REDUCTION OF TRANSIENT DATA TO FREQUENCY RESPONSE 

The transient data measured in flight represent the dy- 
namic response of the airplane to particular control inputs, 
whereas for detailed analysis it is desirable to know the 
airplane response to an arbitrary input. As shown in 
references 1, 2, and 3, it is possible to transform the input 
and output response quantities into frequency-response form. 
This defines the response of the airplane to sinusoidal 
control-surface motion of various frequencies. Since the 
response to any arbitrary input may be obtained by applying 
an inverse transformation to the frequency response, the 


Figure 1 . — Photograph of the test airplane. 
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TABLE II.— FLIGHT-TEST CONDITIONS 


Flight 

number 

Run 

number 

Altitude 

Mach no. 

C L 

W 

c.g. 

KJ 

r 


16 

36, 170 

0. 68 

0. 56 

121, 400 

20. 9 

2. 4 

5. 8 


18a 

35, 500 

. 72 

. 48 

120, 900 

20. 9 

2. 4 

5. 2 


18b 

35, 500 

. 72 

. 48 

120, 900 

20. 9 

2. 4 

5. 2 


20 

36, 260 

. 78 

. 41 

120, 000 

21. 1 

2. 4 

4. 9 

2 

21 

34, 860 

. 82 

. 38 

119, 800 

21. 1 

2. 4 

4. 4 


22 

34, 860 

. 84 

. 33 

119, 300 

21. 1 

2. 4 

4. 3 


23 

34, 100 

. 85 

. 28 

118, 200 

21. 8 

2. 4 

4. 1 


25 

32, 230 

. 63 

. 49 

113, 600 

21. 0 

2. 5 

4. 9 


26 

35, 950 

. 62 

. 61 

113, 200 

21. 0 

2. 5 

5. 8 


15 

35, 070 

. 63 

. 63 

116, 100 

12. 7 

2. 5 

5. 6 


16 

35, 090 

. 66 

. 55 

115, 900 

12. 6 

2. 5 

5. 3 


17 

35, 360 

. 71 

. 48 

115, 900 

12. 6 

2. 5 

5. 0 

o 

18 

35, 060 

. 74 

. 45 

115, 900 

12. 5 

2. 5 

4. 7 

O 

19 

34, 980 

. 77 

. 38 

115, 900 

12. 4 

2. 5 

4. 5 


20 

34, 960 

. 81 

. 37 

115, 700 

12. 5 

2. 5 

4. 4 


21 

35, 150 

. 82 

. 36 

115, 600 

12. 5 

2. 5 

4. 3 


22 

35, 600 

. 84 

. 36 

115, 500 

12. 3 

2. 5 

4. 3 


1 

35, 840 

. 79 

. 40 

119, 100 

24. 9 

2. 4 

4. 8 


3 

34, 170 

. 73 

. 43 

118, 200 

24. 8 

2. 4 

4. 6 


4 

34, 070 

. 69 

. 48 

117, 800 

24. 6 

2. 4 

4. 9 

5 

5 

34, 220 

. 65 

. 54 

117, 000 

24. 6 

2. 4 

5. 3 


6 

35, 100 

. 60 

. 66 

116, 600 

24. 6 

2. 4 

5. 9 


15 

35, 950 

. 72 

. 45 

111, 100 

21. 6 

2. 4 

4. 9 


19 

34, 860 

. 73 

. 43 

116, 200 

29. 7 

2. 4 

4. 8 

A 

22 

34, 590 

. 60 

. 60 

115, 300 

29. 5 

2. 4 

5. 8 

y 

20 

34, 780 

. 70 

. 43 

116, 200 

29. 7 

2. 4 

5. 0 


16 

34, 960 

. 81 

. 37 

116, 800 

29. 8 

2. 4 

4. 4 
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Figure 2. — Two-view drawing of test airplane 


Figure 3. -Typical time histories for elevator pulse and response 

quantities. 
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Figure 4. Typical time histories for aileron pulse and response 
quantities. 
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(a) Initial part of time histories. 

Figure 5. — Typical time histories for rudder pulse and response 
quantities. 
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(b) Complete time histories of p C g B , and /3. 
Figure 5. — Concluded 


frequency response defines all of the basic characteristics 
of the airplane response independently of the particular 
input, and in itself provides a means for studying t lie nature 
of the response and the effect of variables. 

METHOD 

The Fourier integral was used to convert the transient 
responses to frequency-response form. The method used 
for evaluating the Fourier integrals was similar to the one 
described in reference 7 in which ordinates of the time 
history are read at uniform intervals such that the time 
history is closely approximated by a series of parabolic arcs. 
All traces were read at 0.05-second intervals except for the 
pitching velocity, wing deflections, and sideslip angle which 
were read at 0.1-second intervals. The integrals were 
determined by taking an accumulative sum of the product 
of the ordinates and the quadrature coefficients for each 
frequency over the finite length of t lie record. In cases 
where the time history did not end on zero, an analytical 
expression, step or damped sinusoidal oscillation (e. g., 
ref. 8), was used to evaluate the integral from finite time 
to infinity. These calculations were performed on a card- 
programmed IBM 650 digital computer. Corrections were 
made in the data for the dynamic response of the instruments. 

ACCURACY 

The accuracy with which the frequency response may be 
determined is dependent on the relative magnitudes of the 
transform of the transient quantities and the transform of 
errors arising from instrument and reading inaccuracies. 
These errors in the time history were estimated to be of the 
order of 0.005 inch of film deflection. Therefore, for linear 
calibration curves, the transform of the errors consists of the 
Fourier transform of random errors up to about 0.005 inch 
over the length of the transient plus a step error of 0.005 inch 
in the analytical correction at the end of the record. The 
random error was evaluated for a number of traces and was 
found to be generally smaller than the eiror in the end 
correction. For practical purposes the error in the end 
correction may be considered to be a measure of the maxi- 
mum expected error. The data in this report were con- 


sidered to be sufficiently accurate if the amplitude of the 
Fourier transform of a step of 0.005 inch of film deflection 
were less than 10 percent of the amplitude of total transform. 
If the transform of the errors is out of phase with the trans- 
form of the transient, then the error in the amplitude will be 
less than the value above and the phase angle will be in 
error. However, this phase-angle error will be less than 6° 
if the amplitude of the error is less than 10 percent of that of 
the transient. For this reason, accuracy in the amplitude 
also insures accuracy in the phase angles. 

In initial flight tests a number of pilot-applied elevator 
pulses of varying length were recorded and the amplitudes of 
their transforms were compared with the above accuracy 
criterion as shown in figure 6. The error boundary is the 
amplitude of the transform of a step of 0.6° elevator deflec- 
tion (which corresponds to 0.005 inch of film deflection 
multiplied by 10) and represents a boundary of amplitudes 
below which errors of greater than 10 percent can be ex- 
pected. The frequency range was selected to include the 
first two symmetrical structural modes indicated by ground 
vibration tests in reference 9, and the natural frequencies of 
these modes are marked on this figure. It may be seen that 
the transform of the longer pulse (run A) tends to go to zero 
near the frequency of the wing first-bending mode and falls 
below the accuracy boundary at regular intervals thereafter. 
Hence, the longer pulse does not provide adequate excitation 
for accurate evaluation of the frequency response at fre- 
quencies of the structural modes. The shorter pulse (run 
B), on the other hand, provides adequate excitation for fre- 
quencies from 1 to 20 radians per second which is the range 
of interest in this report. Hence, pulses with a time base of 
about 0.5 second, as in run B, were used to obtain the fre- 
quency-response data in this report. A similar procedure 
with aileron and rudder inputs indicated that the primary 
antisymmetric structural modes could be adequately excited 
with pulse inputs. 

REPEATABILITY 

Frequency-response curves evaluated from transient data 
in which the airplane was excited by pulses of different 
lengths at the same flight condition are compared in figure 7. 
Also shown in this figure are boundaries below which errors of 
greater than 10 percent would be expected. This error 
boundary represents the value of d/8 e for which the sum of 
the separate errors of d and 8 e is equal to 10 percent. It may 
be seen that the curves agree within 10 percent over the 
range where the accuracy criterion is satisfied. At fre- 
quencies below 1 radian per second, the frequency response 
obtained by the shorter pulse becomes inaccurate and the 
curves disagree. Also, at frequencies above 7 radians per 
second, the frequency response obtained by the longer pulse 
becomes inaccurate and it may be seen that the values be- 
come erratic. Also shown in this figure is a frequency 
response evaluated from transient data excited by a short 
pulse on another flight in which conditions were slightly 
different. This curve agrees well with the other frequency- 
response curve obtained with a short pulse over the fre- 
quency range of 1 to 20 radians per second, except that it is 
displaced upward at all frequencies. As will be shown later, 
this difference may be accounted for by the difference in the 
time parameter r for the two flights. Thus, the data of 
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Frequency, cj, radians/sec 




Figure 7. — Comparison of pitching-velocity frequency responses 
obtained from short-pulse and long-pulse imputs; M= 0.72. 



LINEARITY 

Since the frequency-response technique is only valid for 
linear systems, frequency responses are questionable when 
obtained at flight conditions in which aerodynamic deriva- 
tives are believed to be nonlinear. The airplane appears to 
have a linear response at 0.72 Mach number and 0.48 lift 
coefficient since the same frequency response was obtained 
with different degrees of excitation in figure 7. However, a 
study of wind-tunnel data indicated that for flight condi- 
tions above 0.8 Mach number and for lift coefficients above 
0.6 at lower Mach numbers the stability derivatives become 
nonlinear. In view of this factor, caution should be exer- 
cised in extrapolating the test results obtained at flight condi- 
tions close to these boundaries (see table II for flight condi- 
tions) to disturbances greater than those used in the flight 
tests, that is, A n max —Y\ g at the center of gravity. 


Figure 6. — Typical elevator control inputs and the amplitudes of their 
Fourier Transforms. 

figure 7 show that with the short-pulse excitation, the 
frequency response can be evaluated with sufficient accuracy 
over the frequency range of interest. 


I. LONGITUDINAL RESPONSE CHARACTERISTICS 

In the following part of the report, the frequency response 
to elevator control-surface motions, as determined from the 
transient flight data, is examined to determine the nature of 
the flexible airplane response. Then the low-frequency part 
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of the response is translated into transfer-function form for 
comparison with values predicted for the rigid and flexible 
airplane. 

MEASURED LONGITUDINAL FREQUENCY RESPONSE 

Longitudinal frequency-response curves evaluated by the 
methods described in Reduction of Transient Data to Fre- 
quency Response are shown in figures 8 and 9 for the test- 
range of Mach numbers at an altitude of 35,000 feet; quanti- 
ties measured at the wing tip and tail as well as at the center 
of gravity are included so that a fairly complete picture of 
the response of the airplane is presented. The frequency 
responses are shown for a forward center-of-gravity location 
in figure 8 and for a rearward one in figure 9. In the follow- 
ing discussion on these frequency-response curves, the various 
longitudinal modes will be identified by comparing the ac- 
celerations at the center of gravity, wing tip, and tail, and 
wing deflections. Also, the frequencies at which these modes 
occur will be compared with those obtained by the ground- 
vibration tests of reference 9. 

'■^riar a 

SHORT-PERIOI) MODE 

The first and largest peak in the amplitude ratio of all of 
the frequency-response curves is the longitudinal short-period 
mode which occurs at frequencies from 1 to 3 radians per 
second. These values agree approximately with frequencies 
estimated from wind-tunnel data. The accelerometer re- 
sponses at wing tip, tail, and center of gravity, which are 
compared in figure 10, show that the accelerations at the three 
locations are essentially in phase, and that the amplitude is 
about the same as the center of gravity and the wing tip, 
but is larger at the tail. Differences between accelerations 
at the various stations are due to contributions of structural 
deflections and pitching acceleration, which will now be 
discussed. 

The wing-bending deflection response, as determined from 
optigraph records, is presented in figure 11. At the short- 
period mode frequency (2 radians/sec), all deflections are 
similarly phased and increase in amplitude toward the wing 
tip, representing a type of deflection similar to the wing 
first-bending mode. Since, for sinusoidal motion, accelera- 
tion is equal to minus the frequency squared times the dis- 
placement, the contribution of the wing deflection to the 
amplitude ratio of the acceleration at the wing tip is of the 
order of —2.5 g per radian. Optigraph records of tail deflec- 
tions indicated that the contribution of tail deflection to the 
amplitude ratio of acceleration at the tail was small at the 
frequency of the short-period mode. 

The contribution of pitching acceleration to the accelera- 
tion amplitude at the wing tip and tail may be simply cal- 
culated by multiplying the amplitude of the pitching 
acceleration (0u>) by the distance to the center of gravity 
(fig. 2). This is permissible because the phase angles of the 
acceleration response at the center of gravity and the pitching 
acceleration response (phase angle of 6/8 equals the phase 
angle of 6/8 plus 90°) are similar at the frequency of 2 radians 
per second. The increment in acceleration response at t lie 
wing tip and tail, then, is about 2.6 and 5.5 g ’ s per radian, 
respectively. Hence, at the wing tip the contributions of 
439029—58 2 


wing deflection and pitching acceleration tend to cancel each 
other so that the acceleration amplitudes at center of gravity 
and wing tip are nearly the same. At the tail, the increase 
in amplitude is primarily caused by the airplane pitching 
acceleration. 

Comparison of the phase angles in figure 11 witli those in 
figure 10 shows that the wing deflections are in phase with 
the accelerations and, therefore, the wing tips deflect in 
proportion to and in the same direction as the airplane 
accelerates. The following sketch illustrates the relationship 
of the wing bending to the airplane center-of-gravity-motion. 
Actually, as shown in the discussion of pitching acceleration, 
the wing tips travel about the same distance in space as the 
center of gravity because of the pitching motion of the 
airplane. 


! 

o 

WING FIRST-BENDING MODE 

The next peak in the frequency response, as seen in figures 
8, 9, 10, and 11, occurs at frequencies of from 8 to 9 radians 
per second, slightly higher than the wing first-bending fre- 
quency (6.9 radians per second) in ground-vibration tests. 
Calculations indicate that the frequency of this mode is higher 
in flight than on the ground because of the additional spring 
force contributed by aerodynamic forces and the increased 
freedom of the body in pitch and translation. As seen in 
figures 8 (a) through (c) and 9 (a) through (c), the response 
peaks for locations on the fuselage are small for this mode as 
compared to those of the short-period mode, but at the wing 
tip (figs. 8 (d) and 9 (d)) a very high peak occurs. Referring 
to figure 10, it is noted that the wing-tip acceleration under- 
goes a 180° phase shift at this peak. Also, according to 
figure 11, the amplitudes of the wing deflections increase to- 
ward the wing tip and are in phase with each other and 180° 
out of phase with the wing-tip acceleration which establishes 
this as the wing first-bending mode. The deflections are also 
in phase with the acceleration at the center of gravity, as was 
the case for the short-period mode. The wing first-bending 
mode is illustrated in the following sketch. The main charac- 
teristic which distinguishes this mode from the short-period 
mode is that the wing deflections are the largest factor in the 
accelerations at the wing tips, while in the short-period mode 
the body translation is the largest factor. 



Several small peaks appear at frequencies from 14 to 16 
radians per second, but these are not well defined because the 
frequency response is in a region of low accuracy. These 
peaks are most prominent on the wing-tip acceleration re- 
sponses (figs. 8 (d) and 9 (d)). The wing-deflection response 
in figure 11 indicates an upward trend in the midsemispan- 
deflection amplitude ratio as compared to that of the wing 
tip. This would indicate a mode of the wing second-bending 
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Frequency, w, radians/sec 


Frequency, cj, radians/sec 


(a) Pitching velocity at the center of gravity. (c) Acceleration at the tail. 

(b) Acceleration at the center of gravity. (d) Acceleration at the w ing tip. 

Figure 8. — frequency responses to elevator inputs for the center of gravity at approximately 12.5-percent c. 
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(a) Pitching velocity at the center of gravity. ( c ) Acceleration at the tail. 

(b) Acceleration at the center of gravity. (d) Acceleration at the wing tip. 

Figure 9.— Frequency responses to elevator inputs for the center of gravity at approximately 2o-percent c. 
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Figure 10. -Comparison of acceleration frequency responses to ele- 
vator input at several stations on the airplane. 


type involving body translation and pitch. Some calcula- 
tions were made on the modes of vibration of the B-47 with 
body translation and pitch included; these indicated a mode 
of the type shown in the following sketch at 16.7 radians per 
second. This agrees qualitatively with the acceleration fre- 



quency responses (figs. 8 (b), (d) and 9 (b), (d)), which show 
that the accelerations of the wing tip and the center of 
gravity tend to become more closely phased when the mode 
becomes prominent at the high Mach numbers. Another 
possible mode of vibration in this frequency range involves 
bending of the inboard nacelle supporting structure which, 
as indicated by ground-vibration tests in reference 9, excites 
considerable wing motion at 16.4 radians per second. This 
mode is probably closely coupled with the mode mentioned 
above, which would explain the appearance of the two closely 
spaced peaks in this frequency region. 

EFFECT OF MACH NUMBER AND DYNAMIC PRESSURE 

Since Mach number and dynamic-pressure effects cannot 
be separated when, as in t he present case, data are available 
for only one altitude, they will be considered together with 



Figure 11. — Wing-deflection frequency responses to elevator input for 
various spanwise stations. 


Mach number arbitrarily selected as the independent 
variable. From figures 8 and 9 it may be seen that there is 
a gradual increase in the amplitude of the short-period mode 
peak for the response quantities up to about a Mach number 
of 0.81, after which the trend reverses, probably because of 
critical Mach number effects. The general level of the 
response quantities also follows this trend at the higher 
frequencies, although the acceleration at the wing tip seems 
to be relatively insensitive to Mach number changes. At 
these higher Mach numbers the frequency response is some- 
what dependent on the magnitude of the elevator pulse 
input because of a nonlinear effect mentioned previously. 

The frequency of the amplitude peak of the short-period 
mode also increases up to a Mach number of 0.81 above 
which it decreases. The frequencies of the peaks of the 
higher modes are relatively constant over this range of 
Mach numbers and dynamic pressures, although the same 
trend may be noted. 

EFFECT OF CENTER-OF-GRAVITY LOCATION 

From figure 12 (a) it is seen that the effect of moving the 
center of gravity from 12.6 to 29.7 lias little effect on the 
frequency response of pitching velocity, although the peaks 
of the short period mode occur at lower frequencies as the 
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(a) Pitching velocity at the center of gravity. 


(b) Acceleration at the center of gravity. 


Figure 12 —Comparison of frequency responses to elevator inputs for various center-of-gravity locations for approximately the same flight 

condition. 


center of gravity is moved back. The effect of center-of- 
gravity movement is more apparent on the acceleration 
response (fig. 12 (b)), which shows a definite trend of the 
short -period mode to higher peaks and lower frequencies as 
the center of gravity is moved back. 

For the rearward center-of-gravity location there are 
from 3,000 to 4,000 pounds of additional fuel in the aft main 
tank (fig. 2) as compared with the forward center-of-gravity 
location. However, there appears to be no significant 
change in the high-frequency portion (5 to 15 radians per 
second) of these responses. This might be expected, since 
the fuselage vertical-bending mode occurs at a much higher 
frequency (29 radians per second). 

EXPERIMENTAL AND PREDICTED TRANSFER FUNCTIONS 

Although frequency-response plots completely defined the 
response characteristics over the frequency range of interest, 
analytical expressions for the response are more useful in 
detailed analyses or in the syntheses of automatic control 
systems. Such analytical expressions, often called transfer 
functions, may be evaluated either from the experimental 


frequency-response plots or from the predicted equations 
of motion of the airplane. In order to show the relationship 
between the transfer function and the frequency response 
more clearly, the predicted transfer functions will be de- 
rived first. Then the method for evaluating experimental 
transfer functions will be explained and, finally, the results 
of the experimental and prediction methods will be compared. 

PREDICTED TRANSFER FUNCTIONS 

Equations of motion. — To define completely the airplane 
dynamic system, the equations of motion must include all 
of the degrees of freedom, but for practical purposes the 
degrees of freedom are usually kept to a minimum by in- 
cluding only the most significant airplane modes. In the 
present analysis, only the short-period mode is considered 
since it was previously shown that the effects of other modes 
on the responses of principal interest, those involving 
motions of the fuselage, were small for the test altitude and 
range of Mach numbers. However, for response quantities 
near the wing tip where at frequencies near and above the 
wing first-bending mode frequency the response is large, 
these equations are, of course, not adequate. 
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The equations for longitudinal motion, where changes in 
lorward velocity are neglected, were used, with pitching- 
acceleration ($) terms added to take account of the distortion 
of the airplane due to pitching-acceleration inertial loads. 
These equations in operator form are: 

[e £a +(C iA +2r)Z»]a+[(r i ,-2r)Z?+^]0+(7 £Se 5 e =O (1) 

( c ma + c V>)«+ [ ( ■„,/)+( c n ,-~ ) z? 2 ] e+c ms 8 „= 0 (2) 

In order to take account of the first-order effects of flexi- 
bility in the preceding equations, the structural deformation 
associated with a coefficient of the variables a or 0 and their 
derivatives was assumed to be in phase with the variable, 
that is, the damping and inertial forces due to structural 
motion were neglected. 

Estimation of stability derivatives. — Stability derivatives 
were derived by available theory with the exception of the 
elevator-effectiveness derivatives, which w^ere modified to 
include additional Mach number effects indicated by wind- 
tunnel tests. The general methods applied were those of 
reference 10 to determine air loads and those of reference 11 
to estimate flexibility effects. Mach number effects were 
included by the Prandtl-Glauert rule as used in reference 10, 
which was indicated by wind-tunnel tests to be not greatly 
in error up to a Mach number of 0.75. 


In applying the general methods mentioned above, several 
modifications were made to facilitate use ol structural test 
data. The structural stiffness of the wing w r as expressed in 
the form of influence coefficients for the front- and rear-spar 
chord wise locations which are compatible with the measured 
influence coefficients and with the spanwise stations used in 
determining the span loading in reference 10. These wing 
influence coefficients were obtained from load deflection data 
presented in reference 12. The weight distribution shown 
in figure 13 wtis separated into equivalent weights as indi- 
cated and used with the influence coefficients for determining 
the wing distortion from inertial loadings. Fuselage influ- 
ence coefficients were obtained from the Boeing Airplane 
Company, the airplane manufacturer. Table III lists all 
the influence coefficients, and Appendix B explains how they 
were used in conjunction with aerodynamic influence co- 
efficients to calculate the lift, moment, and deflection of the 
flexible wing and tail due to any initial angle-of-attack dis- 
tribution. The stability derivatives were calculated by 
determining the initial angle-of-attack distribution due to 
rigid-body motions or to distortion from inertial loads caused 
by rigid-body accelerations, and then calculating the result- 
ing lift and moment coefficients by the method in Appen- 
dix B. The derivatives and the related factors which were 
taken into account are summarized in Appendix C. 

The values of the calculated derivatives are given in table 
IV, which lists separately and aerodynamic and inertial 
contributions. Some of the more important derivatives are 
compared in figure 14 with values calculated for the rigid 
airplane. The Mach number effects are reflected by the 
rigid-airplane derivatives. The flexible-airplane derivatives 
include both Mach number and flexibility effects and show 
that flexibility tends to have the opposite effect of Mach 
number. Of all the derivatives the largest variation occurs 
in C ma , but, in general, aeroelastic effects are not exception- 
ally large because the range of dynamic pressures in these 
tests is not large. 

Evaluation of theoretical transfer functions. — From equa- 
tion (1) and (2), the following transfer function for pitching 
velocity may be obtained: 


8 e 


Kq ( 1 T 7 hD) 


1 + 


2fsp 


D- 


1 


D 2 


and, from the relation n cg =~~ (0—a) 

n cg K n TnJT 2 ) 


1+^ D+^ 2 D 2 

“VsP 


( 3 ) 


( 4 ) 


The acceleration responses at the wing tip and tail are 
obtained by adding the contributions of acceleration at the 
center of gravity, the pitching acceleration, and the struc- 
tural deflection, so that 
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TABLE III.— STRUCTURAL INFLUENCE COEFFICIENTS 

(a) Wing 



1,000-pound load at station 

3Fi 

3 R> 

2 F 

2 R 

IF 

1 R 

5 F 

5 R 


3F 

0. 0659 

0. 0580 

0. 1765 

0. 1663 

0. 2447 

0. 2307 

0. 2033 

0. 1936 


3R 

. 0705 

. 0931 

. 2116 

. 2398 

. 3057 

. 3339 

. 2486 

. 2784 

Inches deflection at 

2F 

. 1769 

. 2052 

. 7419 

. 7795 

1. 2134 

1. 2490 

. 9285 

. 9590 

station 1 

2 It 

. 1804 

. 2318 

. 7881 

. 9005 

1. 3309 

1. 4300 

. 9995 

1. 1070 


IF 

. 2302 

. 2730 

1. 1943 

1. 2934 

2. 3486 

2. 3840 

1. 6374 

1. 7184 


1R 

. 2330 

. 2940 

1. 2430 

1. 4057 

2. 4598 

2. 5920 

1. 7086 

1. 8677 


(b) Fuselage longitudinal coefficients 


Load 

Change in stabilizer angle 

1<7 Normal acceleration _ _ _ 

-0.45° 

( — 0.521° (c.g. = 17 percent), 
\ — 0.506° (c.g. =25 percent) 
+ 0.0342° 

0.000207° 

1 Radian/sec 2 pitching acceleration. _ __ 

1,000-pound tail load (down) _ __ 

1,000 inch-pounds moment applied at stabilizer 


(c) Fuselage lateral-directional coefficients 


Load 

Average change 
in vertical tail 
yaw angle due 
to fuselage and 
tail flexibility, 
deg 

Change in yaw 
angle at tail 
turn meter, deg 

Change in roll 
angle at tail 
turn meter, deg 

1,000 pound aerodynamic load due to side- 
slip of vertical tail 

0. 0685 

0. 0342 

-0. 0798 

1,000 pound aerodynamic load due to rudder 
deflection 

0808 

0402 

. 0942 

1 q lateral acceleration _ 

. 581 

. 290 


1 radian/sec 2 yawing acceleration _ _ 

832 

415 



Note: Wing station locations are shown in figure 13. The letter F denotes front spar location; R denotes rear spar location. 
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TABLE IV.— PREDICTED LONGITUDINAL STABILITY DERIVATIVES FOR THE FLEXIBLE AIRPLANE AT AN 

ALTITUDE OF 35,000 FEET; W= 100,000 POUNDS 


Quantity 

c. g. 


Mach number 



0 

0. 5 

0. 6 

0. 7 

0. 8 



5. 27 

5. 31 

5. 36 

5. 45 

5. 6 

Aerodynamic 

C ^a 


7. l/V 

. 015 

. 013 

. 011 

. 010 

Aerodynamic 

0. 12 

36/ V 

. 069 

. 056 

. 046 

. 038 

° L e 

. 25 

27. 2 /V 

. 050 

. 040 

. 032 

. 025 

Inertial 

. 12 

. 026U 
Q 

. 027 V 
g 

. 348 

. 401 

. 444 

. 473 

C L - e and C K 

. 25 

. 373 

. 431 

. 480 

. 512 

C ^e 

. 12 

-. 034 

-. 031 

-. 030 

029 

028 

. 25 

-. 030 

027 

-. 026 

025 

024 

Cm a 

. 12 

-1. 68 

-1. 81 

- 1. 82 

-1. 82 

-1. 81 

. 25 

-1. 01 

-1. 14 

-1. 15 

-1. 11 

-1. 08 

Aerodynamic 

. 12 

- 27 /r 

053 

-. 047 

-. 040 

- . 036 

^ m a 

. 25 

-26/r 

050 

-. 045 

039 

-. 035 

Aerodynamic 

. 12 

-146/V 

302 

— . 256 

224 

-. 204 

c . 

^ m d 

. 25 

-138 IV 

288 

-. 244 

-. 213 

192 

Inertial 

. 12 

. 024 V 
q 

. 337 

. 398 

. 451 

. 489 

— c m - and c ni 

. 25 

. 024 U 

g 

. 338 

. 400 

. 457 

. 499 

c .. 

. 12 

. 031 

. 030 

. 029 

. 028 

. 028 

L m e 

. 25 

. 030 

. 029 

. 028 

. 027 

. 027 

Cl s 

& e 


. 327 

. 312 

. 306 

. 299 

. 292 

c 

. 12 

-1. 21 

-1. 16 

-1. 14 

-1. 11 

-1. 08 

°e 

. 25 

-1. 17 

-1. 12 

-1. 09 

-1. 07 

-1. 05 

T 


2940 / V 

6. 05 

5. 03 

4. 31 

3. 77 


Since the transfer-function coefficients are complex combi- 
nations of the stability derivatives and mass parameters, 
certain terms have been grouped together so that the coeffi- 
cients are analogous to the familiar equations for rigid 
airplanes. The derivatives C L ^ and C La in the terms C L q 
— 2 t and C L -\- 2 r represent primarily the change in effective 
mass of the airplane resulting from lift due to structural 
distortions. As seen in table IV, the aerodynamic contri- 
butions to these derivatives are relatively small. The 
derivative C m ^ in the term (4tK 2 c/2 V) — C m $ represents the 
change in effective moment of inertia due to pitching mo- 
ments arising from structural distortions induced by rota- 
tional inertial loads. 

Using the stability derivatives given in table IV and 
figure 14 the transfer-function coefficients were calculated 
for both the flexible and rigid airplane and are shown in 
figure 15. The coefficient T n is not shown because its value 

is insignificant. 

Effect of flexibility on transfer-function coefficients. —The 
transfer-function coefficients for the rigid and flexible air- 
plane will now be compared to show for the present tests 
the significance of additional terms included to account for 
flexibility. 


The over-all effect of flexibility on as shown in figure 
15 is to reduce the natural frequency by a maximum of 
about 10 percent. This is principally caused by the reduc- 
tion in C m due to flexibility of the fuselage. This C m 

a « 

effect is partially compensated by a significant increase in 
the —C L C m Q term, which provides an effective spring force 
resulting from the pitching moments arising principally from 
wing deflection due to normal acceleration. The term C m $ 
in the denominator tends to increase the frequency by about 
5 percent. 

The effect of flexibility on damping ratio, is small as 
shown in figure 15. The da ping forces are reduced by 
flexibility, but the spring and i i tial forces are also reduced 
to such an extent that the ov /-all effect of flexibility on 
damping ratio is insignificant. 

The pitching-velocity gai i. ^ of the flexible airplane is 
slightly higher than tnat of t.. igid airplane (fig. 15). This 
change is principally due to tn ' reduction in magnitude of 
C m in the denominator which also is in part compensated 

a 

by the increase in C m ^. Because of this large contribution 
of CmQ for the flexible airplane, increments in gain for a 
10-percent increase in r given in figure 15 show that the 
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Figure 14. — Predicted effect of flexibility on some of the longitudinal 
stability derivatives at an altitude of 35,000 feet. 

flexible airplane has about a 15-percent smaller change in 
gain due to changes in gross weight and altitude. 

Flexibility has little effect on T- as shown in figure 15. 
Although the reduction in C L due to flexibility tends to 
increase T-, this is compensated for the most part by the 
effective reduction in mass in the numerator term. 

The discussion of K- applies to K n , since these factors are 
related by a simple constant. 

The term T n is exceptionally small and is not significant 
over the range of frequencies considered in this report. 

Flexibility tends to increase T n (fig. 15). Although this 
term is small, it does have a noticeable effect at frequencies 
greater than that of the short-period mode. 

Effect of approximate equations on transfer-function coef- 
ficients. — Approximate equations are often used for calcu- 
lating the transfer-function coefficients. These equations, 
which are obtained by simplifying equations 8, 9, 10, and 
11, are as follows: 


— C 






Ke 


~c L a 


La m h. 


Te— 


2 rC m 

2 T 


6V 


(15) 

(16) 

(17) 

(18) 


Transfer-function coefficients for the rigid and flexible 
cases were calculated using both the complete and the ap- 
proximate equations and the results are shown in figure 16. 

Coefficients evaluated by the approximate equations show 
reasonably good agreement with those from the complete 
equations for the rigid airplane, but large discrepancies are 
apparent in the case of u nsp and K- for the flexible airplane. 
This is principally due to the fact that the approximate 
equations reflect only the large reduction in magnitude of 
Cm a due to flexibility while in the complete equations, the 
effect of this reduction in magnitude of Cm a is partially com- 
pensated by a substantial increase in the term —C L C m -. 
Hence, it may be seen that C m ^ should not be neglected in 
calculating the frequency and gain of a flexible airplane. 

EXPERIMENTAL TRANSFER FUNCTIONS 

Experimental transfer functions may be evaluated by 
applying various curve-fitting procedures to the time 
histories or to the frequency responses. Curve fitting of a 
time history may be accomplished by a least-squares fitting 
method such as described in reference 13. Curve fitting of 
the frequency response, as in the follow ing analysis, may be 
done by use of a special set of templates described in refer- 
ence 5. This method may be explained briefly as follows. 
It can be shown that an expression for frequency response 
is obtained by replacing the differential operator, D, in the 
transfer function by the frequency variable, ico. The re- 
sulting complex number can be factored into first- and 
second-order terms expressed in polar (amplitude and phase 
angle) form. The method of reference 5 involves fitting 
the measured frequency response by graphical addition of 
templates selected from a set of curves which represent a 
wide range of first- and second-order factors. 

Evaluation of transfer-function coefficients. — As shown 
previously, the accuracy of the frequency response at low 
frequencies is questionable. For this reason, template 
fitting of the frequency response w^as only considered to be 
valid for frequencies above 1 radian per second. This 
limitation made it difficult to fit the frequency response by 
a unique combination of templates for all of the numerator 
and denominator terms which are involved simultaneously 
in the frequency-response expression. Because of this, the 
natural frequency and damping ratio w r ere determined by 
a least-squares curve-fitting method (similar to that of ref. 
13) of the pitching-velocity time history over the portion of 
the record in which rate of change of elevator position could 
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be neglected. The appropriate templates for these values 
of damping ratio and natural frequency were then fitted to 
the center-of-gravity acceleration frequency-response curves 
to obtain the acceleration gain, K n , and the time constant, 
T nn . The gain K 0 was then determined from K n . By use 
of t his value of K 0 and the previously determined values of 
and a the pitching- velocity frequency response was 
fitted to determine T e . A typical template fit is shown in 
figures 8(a) and (b) for the 0.63 Mach number curve. 

Transfer-function coefficients evaluated from experi- 
mental data for the forward and rearward center-of-gravity 
locations are plotted in figure 17. These values were 
corrected to a common altitude and gross weight in accord- 
ance with predicted variations, but in general these correc- 
tions were small. 

COMPARISON OF EXPERIMENTAL AND PREDICTED TRANSFER FUNCTIONS 

Predicted transfer-function coefficients for the rigid and 
flexible airplane are also shown in figure 17 for comparison 
with the experimental values. The moment of inertia used 
in the predicted coefficients was determined from ground- 
oscillation tests for the basic airplane (ref. 6), with corrections 


made for condition of the fuel tanks. The stability deriva- 
tives of the table IV were employed. 

Effect of Mach number. — Referring to figure 17, it may 
be seen that the experimental values of natural frequency 
show good agreement with the predicted values for the 
flexible airplane and fall about 10 percent below those 
predicted for the rigid airplane. An exception is noted at a 
Mach number of 0.84 where the measured frequency 
decreases sharply, probably due to the large decrease in 
C m as the airplane approaches the pitchup. Also on this 
figure, it is seen that the scatter in the measured damping- 
ratio points prohibits definite confirmation of the predicted 
small variations with Mach number, center-of-gravity 
location, and flexibility. However, the general level of the 
values shows good agreement with theory. 

In figure 17 good agreement is indicated for the pitching- 
velocity gain except at Mach numbers of 0.8 and higher, 
where large unpredicted increases in the measured gain are 
apparent. This again is probably the result of the sudden 
decrease in the static margin as the critical Mach number 
is approached as previously mentioned. On this figure, 
good agreement also is indicated for the pitching-velocity 
time constant, T-, although the experimental values tend 
to fall consistently a small amount below the predicted ones. 


DYNAMIC-RESPONSE CHARACTERISTICS OF A FLEXIBLE 35° SWEPT-WING AIRPLANE 


17 


Effect of center-of-gravity location. — In figure 18, the vari- 
ation of natural frequency and damping ratio with center-of- 
gravity location is shown for several Mach numbers. The 
measured values of natural frequency show good agreement 
with the frequencies predicted for the flexible airplane over 
the test range of center-of-gravity locations. The frequen- 
cies predicted for the rigid airplane are consistently higher 
than those for the flexible airplane, but show about the same 
variation with center-of-gravity movement as for the 
flexible airplane. The measured and predicted damping 
ratios show fairly small changes with center-of-gravity 
location. 

In summary, the longitudinal response of the airplane may 
be adequately predicted by the theory used herein to take 
account of flexibility up to a Mach number of 0.8. The 
largest errors occur in the natural frequency if flexibility is 
not taken into account, but in general all of the transfer- 
function coefficients show only small changes over the air- 
plane range of dynamic pressures at 35,000 feet. 

Since the transfer-function coefficients are relatively 
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Figure 16. — Comparison of longitudinal transfer-function coefficients 
calculated from the approximate equations with those from the com- 
plete equations for the flexible and the rigid airplane; W= 100,000 
pounds, c. g. ~= 25-percent M. A. C. 



M 

Figure 17. —Comparison of predicted and experimental values of 
longitudinal transfer-function coefficients; 1^=115,000 pounds. 


insensitive to flexibility effects at these flight conditions, a 
more sensitive transfer function, the wing-tip deflection 
response, was calculated from equation (7) in order to pro- 
vide a more accurate check on the aeroelastic calculations. 
The method of calculating deflections is shown in Appendix 
B. The comparison between experimental and predicted 


18 


REPORT 1330 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Center of gravity, percent M.A.C. 


Figure 18.— Comparison of predicted and experimental variation with 
center-of-gravity position of natural frequency and damping ratio of 
the short-period mode; W= 115,000 pounds. 

values in frequency-response form is shown in figure 19. It 
may be seen that the amplitude ratios agree very well and 
that the phase angles of the experimental response tend to 
lag the predicted ones 20° to 30° up to frequencies of 5 
radians per second. From this, it would appear that the 
theoretical calculations of the contribution of the wing to 
the stability derivatives are accurate except for a small 
phase lag which is probably due to the inertial and damping 
forces due to wing bending motion which were neglected in 
the theo^. 

II. LATERAL-DIRECTIONAL RESPONSE 
CHARACTERISTICS 

In the following part of the report, the equations of 
lateral-directional motion are presented first to give the 
reader the background needed to interpret dynamic response 
characteristics in terms of transfer functions and stability 
derivatives. Frequency responses determined from tran- 
sient flight data are then examined with this knowledge and 
comparisons are made between predicted flexible and rigid 
airplane transfer -function coefficients and those evaluated 
from flight data. Finally, aerodynamic derivatives extracted 
from the transfer-function coefficients are compared with 
es t im ated d eriv a t i ves . 

PREDICTED RESPONSES 

As in the longitudinal case, considerable simplification in 
the equations of motion can be realized if the structural 
deformations are assumed to be in phase with the aero- 
dynamic and inertial loads resulting from airplane motions. 



Figure 19. — Comparison of experimental and predicted wing-tip-de 
flection frequency response to elevator input; flight 3, run 17. 


The aerodynamic forces arising from these deformations are 
included by modifying the derivatives in t lie lateral-direc- 
tional equations of motion for a rigid airplane. 

EQUATIONS OF MOTION 

Three-degree-of-freedom lateral-directional equations of 
motion with respect to stability axes are given in Appendix 
D. For the flexible airplane appropriate derivatives cor- 
responding to structural deflections due to inertial and dead- 
weight forces resulting from airplane motions have been 
added to the equations. These additional terms are then 
combined with the corresponding rigid airplane inertias to 
form an effective inertia for the flexible airplane. As 
explained in Appendix D, it was found to be permissible to 
neglect a number of the inertial derivatives since the cumula- 
tive effect of these quantities on the transfer-function 
coefficients was found to be negligible. 

AERODYNAMIC DERIVATIVES 

The stability derivatives were obtained by available 
theory and wind-tunnel data. Effects of structural deflec- 
tions on the derivatives were obtained through use of 
aerodynamic and structural influence coefficients in the 


DYNAMIC-RESPONSE CHARACTERISTICS OF A FLEXIBLE 35° SWEPT-WING AIRPLANE 


19 


maimer described in Appendix B. The chordwise centers of 
pressure of all aerodynamic loads, except loads due to the 
rudder and ailerons, were assumed to be at the quarter-chord 
line of their respective surfaces. The average chordwise 
center of pressure for loads due to the rudder was estimated 
to be the 0.48 chord from the pressure-distribution data given 
in reference 14. Although no pressure-distribution data 
were available for the specific wing-aileron plan form, the 
center of pressure for the aileron loads was estimated to be 
at t lie 0.58 chord from an examination of pressure distribu- 


tions on similar swept-wing plan forms. Fuselage influence 
coefficients used in the analysis were obtained from data 
given in references 15 and 16 and are listed in table III 
together with the wing influence coefficients which were 
obtained from reference 12. Wing stations and equivalent 
weights used in the analysis are shown in figure 13. An 
example of the application of the method to obtain the 
effect of wing deflections on the derivative is given in 
Appendix B. Predicted aerodynamic derivatives for both 
the flexible and rigid conditions are listed in table V. 


TABLE V— PREDICTED LATERAL-DIRECTIONAL STABILITY DERIVATIVES FOR THE FLEXIBLE AND RIGID AIRPLANE 

AT AN ALTITUDE OF 35,000 FEET; 1V= 115,000 POUNDS 


Quantity 

Flexible 

Rigid 

M 

0. 5 

0. 6 

0. 7 

0. 8 

0. 5 

0. 6 

0. 7 

0. 8 

Q 

87. 2 

125. 5 

171. 0 

223. 1 

87. 2 

125. 5 

171. 0 

223. 1 

C L 

. 922 

. 642 

. 471 

. 361 

. 922 

. 642 

. 471 

. 361 


-. 430 

-. 429 

-. 427 

-. 426 

-. 490 

507 

-. 540 

-. 570 

c„ 

p 

-. 0960 

0660 

-. 0470 

-. 0350 

-. 0994 

0692 

-. 0507 

-. 0389 

C>r 

. 215 

. 136 

. 090 

. 060 

. 264 

. 184 

. 135 

. 104 

Cn r 

-. 142 

137 

-. 134 

-. 130 

150 

-. 150 

-. 150 

-. 150 

C 'e 

-. 170 

-. 141 

-. 122 

-. 112 

-. 157 

129 

113 

-. 101 

Cn e 

. 110 

. 106 

. 102 

. 098 

. 120 

. 120 

. 120 

. 120 

Cy <3 

-. 575 

-. 566 

-. 555 

-. 544 

600 

600 

-. 600 

-. 600 

C \ 

. 0052 

. 0090 

. 0109 

. 0122 

. 0057 

. 0095 

. 0119 

. 0135 


-. 0724 

0700 

0672 

-. 0645 

-. 0778 

-. 0778 

-. 0778 

-. 0778 


. 183 

. 180 

. 175 

. 171 

. 194 

. 194 

. 194 

. 194 

c '<. 

. 0594 

. 0528 

. 0462 

. 0392 

. 0770 

. 0770 

. 0770 

. 0770 


-. 0103 

0064 

0041 

0027 

-. 0133 

0093 

-. 0068 

-. 0052 

r* 

-. 116 

059 

005 

. 031 

-. 103 

047 

-. 010 

. 014 

rz 

-. 048 

021 

005 

. 007 

047 

020 

-. 005 

. 007 


TRANSFER-FUNCTION COEFFICIENTS 

Transfer functions obtained by solving the three-degree- 
of-freedom equations are presented in Appendix D. As 
noted in this appendix, some of the transfer functions may 
be simplified for the frequency range of interest. After the 
spiral mode is neglected, rolling response to aileron becomes 

p_ Az(D 2 -\-diD-\-a2) 

8 a ==s (D+D r )(D 2 +c l D+c 2 ) 


d T D 2 +CiD+C 2 

These quantities may also be expressed as 

r K r/ s r D 

8 r l-\-2t D liD/u nDR -\-(D/Un DR ) 2 



8 r l+2£ DR D/u nDR +(D/u nDR ) 2 


(22) 


(23) 

(24) 


The transfer function may also be written as 

P Kp, * fl [l + 2 £ a D/ 03 a + (Z)/o> a ) 2 ] 

8 a (1 + 77)) [l+2£ DR D/a) nDR +(D/o) nDR y] 


(20) 


For yawing and sideslipping response to rudder motion, both 
spiral and rolling modes may be neglected so that 


H r D 

8r D 2 + CiD + C 2 


( 21 ) 


The predicted frequency responses may be obtained from 
the transfer functions after io: is substituted for the operator 
D. The conversion of numerator transfer-function coeffi- 
cients from stability to body axes for comparison with fre- 
quency responses measured with respect to body axes is 
given in Appendix D. 

APPROXIMATE EQUATIONS TO OBTAIN AERODYNAMIC DERIVATIVES 
FROM TRANSFER-FUNCTION COEFFICIENTS 

The equations for the aerodynamic derivatives were ob- 
tained as follows: The coefficients C z , C 2) and Ci of the char- 
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acteristic equation given in Appendix D were expressed in 
terms of their principal aerodynamic derivatives. The 
coefficients were also expressed in terms of the factors of the 
characteristic equation D r , c u and c 2 . By equating the 
coefficients of like powers of D, three simultaneous equations 
were obtained with three unknowns, L p , N r -\-Yp, and Np. 
A cubic equation for L p was then obtained by eliminating 
the other unknowns. The derivative L p was found by using 
an approximate solution for a cubic equation together with 
an approximate equation for the amplitude of the roll to 


Equations for the 


yaw^ ratio of the Dutch roll mode, 

two other derivatives could then be obtained. The equations 
are 


L p — T) r n 
N r -\- Yp=fx — Ci 
Ne'=c 2 +»(Dr- Cl ) 


where 



(N p —Y p —D r r z ) 

Vl.+D r 2 /c 2 


(25) 

(26) 
(27) 


(28) 


Note that these equations correspond to those given in refer- 
ence 3 except for the quantity /z- The additional term /z 
could be considered as a measure of the coupling between the 
Dutch roll and rolling modes. The equation used for Lp 
was obtained from reference 3 


L d ' = 


Ntj'\l-\-L p 2 /c2 


(29) 


Approximate relations for control effectiveness were obtained 
by substituting typical numerical values in the equations 
given in Appendix D, and the following expressions (anal- 
ogous to those given in ref. 3) w r ere found to be valid for the 
present airplane. The rudder effectiveness derivative w 7 as 
obtained from 


N b =H r or — Jp (30) 

while aileron effectiveness w r as obtained from 







M 


Figure 20. Influence of flexibility on effective lateral-directional 
inertia characteristics. 


function coefficients were first calculated from the predicted 
derivatives (table V) using the complete expressions for the 
transfer-function coefficients. Then approximate values of 
the derivatives were calculated from the predicted transfer- 
function coefficients using equations (25) to (29), and also 
using the equations with ju=0, which correspond to the equa- 
tions given in reference 3. Comparisons of these deriva- 
tives are shown in figure 21. Also shown are derivatives 
calculated from equations given in reference 3. In all cases 
the values of derivatives using the approximate equations 
are closer to the estimated values than those obtained using 
the equations from reference 3. This is particularly true for 
the derivative, C nr +2K z 2 C Y p- Although not shown, a 
similar comparison was made for the rigid case, and the 
inclusion of the quantity /x also improved the agreement with 
the predicted derivatives. 




(31) 


MEASURED RESPONSES AND COMPARISON WITH 
PREDICTED RESPONSES 


Aerodynamic derivatives may be found from the quantities 
on the left-hand side of equations (25) through (31) through 
use of the definitions given in the notation. 

Thus the aerodynamic derivatives may be obtained from 
transfer-function coefficients matched from experimental 
data, provided that an estimate is made of the product of 
inertia, the derivative C n , and also the effective moments of 
inertia for the flexible airplane. Estimated variations in 
effective inertia characteristics are presented in figure 20. 
The principal effect of flexibility is to reduce the rolling 
moment of inertia. This effect is largely due to the mass of 
the wing-mounted nacelles. In order to demonstrate the 
adequacy of the approximate equations, predicted transfer- 


FREQUENCY RESPONSES 

In the discussion that follows, effects of the relatively low 
frequency lateral-directional modes will be investigated first. 
Comparisons will be made between measured and predicted 
responses to determine the frequency range for which the 
predictions are adequate. Then the higher frequency 
structural modes will be discussed and compared with values 
obtained from the ground oscillation tests of reference 9. 
The frequenc 3 ^-response data to be presented include: 
measured frequency responses to aileron and rudder inputs 
over a fairly wide Mach number range (figs. 22 and 23), 
typical comparisons of measured responses at the center of 
gravity and tail (figs. 24 and 25), and typical comparisons 
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of measured and estimated frequency responses (figs. 26 
and 27). 

Lateral-directional modes. — The ailerons excite principally 
the rolling mode as well as some Dutch roll mode as shown 
in the roll-rate response (fig. 22 (a)), which is seen to be 
of the same form as that indicated in equation (20). The 
rolling mode (a first-order term) causes the net shift of 90° in 
the phase angle between the frequencies of 0.5 to 5 radians 
per second and the slope of the amplitude curve in the 
frequency range of 2 to 7 radians per second, while the peak 
in the amplitude curve that occurs at frequencies from about 
1 to 1.5 radians per second is due to the Dutch roll mode 
(a second-order term). Responses in sideslip and yaw rate 
are quite small (figs. 22 (b) and 22 (c)). The rudder mainly 
excites the Dutch roll mode (fig. 23). For both inputs, 
differences in tail and center-of-gravity roll and yaw rates 
are small in this frequency range. This may also be seen from 
the time histories (figs. 4 and 5). 

The predicted results for the flexible airplane in general 
agree reasonably well with experimental results for fre- 
quencies up to about 3 to 8 radians per second, with the 
possible exception of wing-tip acceleration (figs. 26 (e) and 
27 (f)). Also some discrepancies are apparent in the lowest 
frequency range (near 0.5 radian/sec) which is the region in 
which the experimental frequency -response accuracy ob- 
tained from pulse-type inputs becomes relatively poor, as 
was discussed in the section Reduction of Transient Data to 
Frequency Response. For the aileron responses, the pre- 
dicted gain for the flexible case is much lower than that for 
the rigid, reflecting the reduction in aileron effectiveness. 
However, the p/5 r responses (fig. 27 (a)) indicate a higher 
gain for the flexible case. This is largely due to the higher 
value of the ratio of C i0 to the effective moment of inertia in 
the flexible case. Flexibility also reduces the frequency and 
damping of the Dutch roll mode. A discrepancy between 
theory and experiment in the curves for r/8 a (fig. 26 (b)) is 
that the experimental phase angles shift 180° while the pre- 
dicted values approach a 360° shift. One possible explana- 
tion is that the smaller phase shift could have been obtained 
from the predicted curve if a more negative value of C n6a had 

been assumed. A structural mode could also affect the 
phase angle in this region. 

Structural modes. — A peak occurs in all the frequency 
responses to aileron (fig. 22) at about 13.4 radians per second 
which will be called the first wing antisymmetric bending 
mode. This compares with 14.4 radians per second obtained 
from ground vibration tests (ref. 9). The mode is partic- 
ularly noticeable in the center-of-gravity roll rate and wing- 
tip acceleration responses. In the time history for the wing 
optigraph, this mode is seen to predominate. Another point 
of interest is that the amplitude of the yaw rate at the tail is 
greater than that for the center of gravity, as may be seen 
from both the frequency-response and time-history curves. 
The mode frequency remains essentially constant for the 
flight range investigated although there is a slight increase in 
amplitude with Mach number. Measurements of the time 
history indicated that the damping was also essentially 
constant over the flight range at a value of about 0.025 
critical. 





Derivatives calculated by equations from this report. 

Derivatives calculated by equations from reference 3. 



M 

Figure 21. — Comparison of predicted lateral-directional stability de- 
rivatives for the flexible airplane with those calculated from predicted 
transfer-function coefficients using approximate equations. 


The shape of this mode, obtained from measurements 
of the free oscillation amplitudes of the accelerometer and 
turn-meter traces resulting from aileron inputs, is shown in 
figure 28. When the curve was determined, small correc- 
tions for chordwise displacement of the instruments were 
made where necessary to transfer the accelerometer results 
to the 17-percent and 58-percent chord lines corresponding 
to the front and rear spars, respectively. In the placement 
of pickups for an autopilot system, it is desirable to locate 
them at points on the structure that are least affected by 
the structural modes in order to avoid introducing undesired 
signals. For a roll gyro, the optimum location for minimiz- 
ing the angular velocity input of the antisymmetric wing 
mode would be where the wing mode shape has zero slope 
which, as shown in figure 28, is at rj=0.3S near the inboard 
nacelle, while for a linear accelerometer, the optimum loca- 
tion would be at a nodal point, which in figure 28 is at rj of 
about 0.55 just outboard of the nacelle. The wing acceler- 
ometer node line obtained from the two span positions 
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ter of gravity. center of gravity. 


Figure 22. — Frequency responses to aileron inputs at various Mach numbers. 




Frequency, w, radians/sec 


(a) Rolling velocity at the center of gravity. (b) Yawing velocity at the center of gravity. (c) Yawing velocity at the tail. 

(d) Sideslip angle. (e) Wing-tip acceleration. 


Figure 23. — Frequency responses to rudder inputs at various Mach numbers. 
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Figure 24. — Comparison of the measured yawing velocity frequency 
responses at the center of gravity and the tail to an aileron input; 
Af=0.71. 


shows good agreement with that obtained from ground 
vibration tests. 

In the wing-tip acceleration curve (fig. 22 (d)), a dip in 
the amplitude and a shift of approximately 180° in the 
phase angle relative to the center-of-gravity roll rate (fig. 
22 (a)) occurs in a frequency range somewhat lower than 
the wing bending frequency. This motion can be illus- 
trated as follows: At low frequencies, the relation between 
wing bending and roll angle is as shown in the following 
sketch. As the frequency increases, the amount of wing 



bending relative to roll angle becomes greater so that the 
acceleration at the tip reverses sign with respect to the roll 
angle as shown in the next sketch. This trend continues 



until the bending frequency is reached. 

An examination of figure 22 shows that no additional 
modes were excited to any appreciable extent in the measured 
frequency range above that of the wing first antisymmetric 
bending mode, although two modes were measured from 
ground vibration tests (ref. 9). 

While the rudder pulses also excite the wing anti- 
symmetric bending mode, their principal effect is to excite 
two slightly higher frequency modes which are predominantly 



(a) Rolling velocity. 

(b) Yawing velocity. 

Figure 25.— Comparison of measured frequency responses at the 
center of gravity and the tail to a rudder input; Af=0.7l. 

fuselage modes as may be seen from the responses of yaw 
and roll rate at the tail and center of gravity (figs. 23 and 
25). The lower frequency, which is predominantly fuselage 
side bending, is 16.5 radians per second while the upper 
frequency, which is principally fuselage torsion, is 23 radians 
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(a) Rolling velocity at the center of gravity. (b) Yawing velocity at the center of gravity. (c) Yawing velocity at the tail. 

(d) Sideslip angle. (e) Wing-tip acceleration. 

Figure 26. — Comparison of experimental and predicted frequency responses to aileron input; M-- 0.71. 
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(a) Rolling velocity at the center of gravity. (b) Rolling velocity at the tail. (c) Yawing velocity at the center of gravity, 

(d) Yawing velocity at the tail. (e) Sideslip angle. (f) Wing-tip acceleration. 


Figure 27. — Comparison of experimental and predicted frequency responses to rudder input; M= 0.71. 
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Figure 28. — First antisymmetric wing bending free oscillation mode 
shape; M= 0.71. 


per second. These frequencies compare with ground vibra- 
tion test values of 18.0 and 24.5 radians per second, respec- 
tively (ref. 9). The prominence of these modes at the tail 
may also be seen from the transient responses due to the 
rudder input (fig. 5). 

The frequency responses for roll and yaw rate at the tail 
due to a rudder pulse (fig. 25) also indicate regions of low 
amplitude at frequencies below the structural natural 
frequencies (about 3 to 8 radians per second). While phase 
angle for the tail yaw rate shifts 180° out of phase with that 
of the center of gravity, as was the case for the wing at 
frequencies below that of the antisymmetric bending mode, 
the phase angle shifts for the tail and center-of-gravity roll 
rate are approximately the same so that they remain in 
phase. It is of interest to note that by the inclusion of only 
the pseudostatic effects of sideslip, yawing acceleration, etc., 
on the tail deflections, these amplitude and phase-angle 
trends are predicted (figs. 27 (b) and 27 (d)). While the 
direction of the 180° phase-angle shift does not agree for the 
tail yaw rate, the resultant shift is the same. Since the 
amplitude ratio is quite low in this region, some question 
does exist as to the direction of the phase-angle shifts for 
both the experimental and predicted values. 

TRANSFER-FUNCTION COEFFICIENTS 

Transfer-function coefficients of the simplified transfer 
functions for p/8 a , r/d r , and 0/5 r described in a previous section 
were obtained by matching the response time histories (e. g., 
figs. 4 and 5) by means of an analog computer. They were 
then converted from body axes to stability axes. In some 
cases, small residual oscillations, which were principally due 
to the lightly damped Dutch roll mode, were present when 
the airplane control pulse was applied. Initial conditions 
were applied to the simulator to include this motion, although 
in all cases this had only a small effect on transfer-function 
coefficients. In matching the time histories of p/8 a on the 
analog computer, it was found convenient to represent the 
third-order transfer function by parallel networks of first- 
order and second-order transfer functions. This representa- 
tion is equivalent to separating the third-order transfer 
function into the partial fraction form shown in Appendix D. 
In this way contributions of the Dutch roll and rolling modes 
to the resultant motions could be matched separately. 


Comparisons of experimentally determined transfer- 
function coefficients with predicted values (with respect to 
stability axes) for both the flexible and rigid case based on 
the simplified transfer functions (eqs. (20), (23), and (24)) 
are shown in figure 29. In general, agreement between 
experiment and theory lor both numerator and denominator 
is seen to be good, as was the case reported in reference 17 
for 25,000 feet. 

For the denominator terms, the values of measured damp- 
ing ratio are an average of 0.015 less than the estimated 
flexible values. While the ratios of measured to estimated 
values are relatively large, the average magnitude is ap- 
proaching zero at the lower values of Mach number, and 
hence the discrepancy is considered to be small. Also some 
scatter occurred in the time constant T. This parameter was 
difficult to evaluate because it did not have as important an ef- 
fect on the transient response to a pulse as did the other param- 
eters. Flexibility is seen to reduce both the frequency and 
damping of the Dutch roll mode. These reductions are 
caused primarily by the decrease in vertical-tail effectiveness 
due to luselage and tail flexibility having a larger effect than 
the decrease in effective yawing moment of inertia. Flexi- 
bility changes the time constant T very little. This term 
depends principally on the ratio of effective rolling moment 
of inertia to the aerodynamic damping in roll. Since the 
rolling moment of inertia is due largely to the nacelle and 
wing weights, wing flexibility reduces both the aerodynamic 
and inertia loads by about the same amount. 

For the numerator terms also, good agreement is obtained 
between theory and experiment. The scatter that occurred 
in the measured values of resulted from the small effect 
of fa on the time histories matched with the analog computer. 
The comparison of the matched curves in figures 4 and 5 
with the measured time histories indicates that the transfer 
functions corresponding to the Dutch roll mode only tor the 
rudder pulses, and the rolling and Dutch roll modes for the 
aileron pulses were sufficient to define the time histories ade- 
quately. The greatest effect of flexibility is seen to be in 
the reduction in gain for p/8 a . This reflects the reduction 
in aileron effectiveness associated with outboard ailerons. 
The predicted values of and co a are slightly different from 
$ DR and u nDR and approach these quantities at the higher 
M ach numbers. These differences depend upon the aero- 
dynamic derivatives that affect the quantities /x and C n ^, 

which become relatively small at the higher Mach numbers 
and corresponding low angles of attack. These effects are 
reflected in the time histories by the amount of excitation 
of the Dutch roll mode and also by a reduction in over-all 
gain. For the limiting case of { a =t DR and a> a =co„ D/2 , effects 
of the Dutch roll mode are eliminated and the transfer 
function (eq. (20)) reduces to a first-order form. 

AERODYNAMIC DERIVATIVES 

Generally good agreement is obtained between predicted 
and experimental values of the aerodynamic derivatives 
(fig. 30). Differences between predicted aerodynamic de- 
rivatives for the flexible and rigid cases are generally some- 
what greater than the differences between transfer-function 
coefficients since the latter are partially compensated by 
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Figure 29. Comparison of experimental and predicted lateral-direc- 
tional transfer-function coefficients. 
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Figure 30. — Comparison of experimental and predicted lateral-direc- 
tional stability derivatives. 


DYNAMIC-RESPONSE CHARACTERISTICS OF A FLEXIBLE 35° SWEPT-WING AIRPLANE 


29 


changes in the effective moments of inertia due to flexibility. 
In general, flexibility tends to reduce the values of the aero- 
dynamic derivatives with the largest reduction occurring 
for C,, and C, . However, flexibility tends to increase the 

magnitude of C tf} since the usual decrease in value of a 
derivative due to flexibility is more than compensated by 
the increase in dihedral due to wing lift. 

CONCLUSIONS 

The evaluation of the dynamic response of a large flexible 
airplane to elevator, aileron, and rudder pulses at 35,000 
feet over a Mach number range of 0.6 to 0.85 and the com- 
parisons with predicted response have led to the following 
conclusions : 

1. The pulse technique provided dynamic-response data 
which were sufficiently accurate for evaluation of frequency 
response from 1 to 25 radians per second. This frequency 
range includes the relatively low frequency longitudinal 
short-period, Dutch roll, and rolling modes, as well as the 
higher frequency structural modes, wing first symmetric 
bending, wing first antisymmetric bending, fuselage vertical 
bending, fuselage side bending, and fuselage torsion modes. 

2. Predicted longitudinal and lateral-directional responses 
based on rigid airplane equations witli coefficients modified 
to include zero-frequency aeroelastic effects were in good 
agreement with measured responses for frequencies below 
one-half of the lowest structural mode natural frequency. 

3. For longitudinal responses, the principal effect of air- 
plane flexibility is to decrease the natural frequency of the 
short -period mode, a trend which may be attributed prin- 
cipally to fuselage and tail bending and the associated loss 
in the angle-of-attack moment derivative, C ma . This C ma 


effect would be more significant and have a larger effect on 
the steady maneuvering acceleration per unit elevator de- 
flection if it were not for the compensating effect of deflec- 
tions due to inertial loads. 

4. The principal effect of flexibility on the lateral-direc- 
tional response is the reduction in gain of airplane response 
to aileron motion, reflecting the reduction in aileron effective- 
ness, and also a reduction in damping and frequency of the 
Dutch roll mode. 

5. The principal lateral-directional and control-effective- 
ness derivatives evaluated with approximate equations were 
in good agreement with predicted derivatives when zero- 
frequency effects of flexibility were included. 

6. For frequencies near and above the natural frequencies 
of structural modes, the modified equations for zero-fre- 
quency effects of flexibility were inadequate. However, 
near the natural frequency of a structural mode, the re- 
sponse amplitude varies considerably with location on the 
airplane, and it appears that for some locations, at least, the 
response near the structural mode frequency would be the 
same as predicted by the simplified equations. 

7. The frequencies of the structural modes measured in 
flight were within about 10 percent of the frequencies of the 
corresponding modes measured in ground vibration tests. 

8. The method of employing aerodynamic and structural 
influence coefficients in the aeroelastic calculations proved 
to be advantageous in simplifying spanwise loading calcula- 
tions and in applying static test load data. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Autj. 6 , 1957 
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NOTATION 


aspect ratio 
lift coefficient 
rolling-moment coefficient 
pitching-moment coefficient 
yawing-moment coefficient 
side-force coefficient 


Ci 

9 

i 

m 

m F 


differential operator. 

at 

moment of inertia about the axis, 
slug-ft 2 

effective rolling moment of inertia for 
flexible airplane, I x —0 o SbCi », slug-ft 2 
moment of inertia about the 7 axis, 
slug-ft 2 

effective yawing moment of inertia for 
flexible airplane, I z —ci 0 SbC tl - , slug-ft 2 
product of inertia, slug-ft 2 
dimensionless radius of gyration about 
principal lateral axis, mean aerody- 
namic chords 

dimensionless radius of gyration about 
principal normal axis, wing spans 


n 

V 

<Z 

Qo 

r 

t 

x { ) 


V 

z 


71 

gain of — transfer function 

o e 


gain of simplified 7- transfer function 

Oa 


Zo 

a 

P 


gain of simplified — transfer function 
0 T 


8a 


gain of simplified — transfer function 
1 8 r 


gain of — transfer function 
8 e 


8e 

8 r 


Mach number 

wing area, sq ft 

rolling-mode time constant, sec 

normal acceleration time constant, sec 

pitching velocity time constant, sec 

velocity, ft/sec 

gross weight of airplane, lb 

wing span, ft 

local chord of wing measured paraded to 
the plane of symmetry 
center of gravity, percent c 
wing mean aerodynamic chord, 


f 

fa 

V 

6 

M 



section lift coefficient 

acceleration due to gravity, 32.2 ft/sec 2 

v=i 

mass of airplane, slugs 

effective mass for lateral acceleration of 


flexible airplane, m—- 0 ^ C y . 

V 8 

normal acceleration, positive downward, 
gravity units 

rolling velocity, radians/sec 
dynamic pressure, lb/sq in. 
dynamic pressure, lb/sq ft 
yawing velocity, radians/sec 
time, sec 

longitudinal distance from center of 
gravity to subscript quantity, positive 
when center of gravity is forward of 
subscript quantity, ft 
spanwise coordinate perpendicular to 
plane of symmetry, ft 
structural deflection, positive downward, 
in. 

vertical coordinate, positive downward, 
in. 

angle of attack, radians 
angle of sideslip, radians (except as other- 
wise noted) 

total aileron deflection, measured in a 
plane perpendicular to the hinge line, 
positive right aileron up, radians (ex- 
cept as otherwise noted) 
elevator deflection, positive downward, 
radians (except as otherwise noted) 
rudder deflection, measured in a plane 
perpendicular to the hinge line, posi- 
tive trailing edge left, radians (except 
as otherwise noted) 

change in streamwise angle of attack due 
to wing distortion, radians 
damping ratio, dimensionless 

damping ratio of numerator term of ~ 

8a 

transfer function 

spanwise coordinate, wing semispans 
pitch angle, radians 

parameter used in evaluating aerody- 
namic derivatives from, transfer- func- 
tion coefficients, defined in equation 
( 28 ) 
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mass density of air, slugs/cu ft 

m 

Jsv ’ scc 
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/ output \ 
\ input / 


CO 

co„ 


angle of bank, radians 
phase angle of output quantity minus 
phase angle of input quantity 
angle of yaw, radians 
frequency, radians/sec 
undamped natural frequency of numer- 
ator term of ^ transfer function 
8 a 

undamped natural frequency, radians/sec 


c La , c L& , c l . 
S c ma 

Cm u , C m ^ } 
r r r 

n rt p 

^ n g 

S- S’ S 

S. s 

c,. t Ci.., r 

<p t + 


derivative of coefficient with respect to 
> subscript 


Ci r , 0 nr , Ci p , Cn p derivative of coefficient with respect to 
subscript Xyp 


g oSb 2 r 

OT 7T Ll n’ 


2 VIx 


per sec 


L, 


q 0 Sb 2 n 

2VU° h ’ PPr 


q 0 Sb n 2 

V— Cl , per sec 2 

1 x F 


N 

1 v n 


q 0 Sb 2 

2 Viz, 


C n , per sec 


N r 

q 0 Sb* n 
2 VI , r °- ^ 

Nf> 

(] 0 & b x-y 2 

j C n per sec 2 
±z F 

Y, 

QoS n 

ir c v>' crsec 

Y, 

g 

per sec 

Lb 

(1 qS b 2 

— Ci v per sec 2 

1 * F 

Ns 

q„Sb n 2 

r Cn V per sec 2 
lz F 

Y s 

q 0 S 
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IxzJ-° Sh0 'i 
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Lp+rxpNp 
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Lp-\-rx F Np 

N/ 

N r -\-r ZF L r 

Nff' 

Nfi-\-r Zp Lp 

u 

Ld-^rxpNs 

Ns’ 

Nb~\r r z F Ls 


SUBSCRIPTS 

B 

body axis 

DR 

Dutch roll 

eg 

center of gravity 

SP 

short period 

t 

tail 

wt 

left wing tip 


Dots are used to indicate differentiation with respect to 
da 


time; for example, a= 


dt 
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A METHOD FOR DETERMINING THE AERODYNAMIC LIFT 
AND MOMENT OF A FLEXIBLE WING THROUGH USE OF 
INFLUENCE COEFFICIENTS 

The approach used is generally that presented in reference 
11 . However, for the calculation of wing deformations, the 
use of structural influence coefficients in conjunction with 
aerodynamic influence coefficients will be introduced in 
place of the wing bending and torsional stiffness distributions 
El and GJ, and the distributed aerodynamic loadings used 
in reference 11 . 

LONGITUDINAL STABILITY DERIVATIVES 

The loading on a flexible wing may be separated into two 
parts: that of the rigid wing and that produced bv the 
wing deflection. The discussion that follows will be con- 
cerned with the determination of the loading produced by 
a wing deflection corresponding to a given rigid -wing 
loading. If a relaxation procedure is used in which suc- 
cessive deflections are calculated from each aerodynamic 
load (i. e., the first deflection due to the rigid wing aero- 
dynamic loading, the second deflection due to the aerody- 
namic loading resulting from the first deflection and so on), 
then the change in streamwise angle of attack may be 
written as follows: 


e (v) — €oW#+€i(^)<Z 1 2 +€ 2 ( 77 )# 3 + ... (Bl) 

where 

€ 0 ( 77 ) is the angle-of-attack distribution due to rigid-wing 
loading 

61(77) is the angle-of-attack distribution due to loading 

obtained from 60 ( 77 ) 

62 ( 77 ) is the angle-of-attack distribution due to loading 

obtained from 6 ^ 77 ) 

The incremental angle-of-attack distributions 60 ( 77 ), 61 ( 77 ) 

may be converted into incremental lifts to form the series 

*C L =C LA q+C LB q*+C Lc q*+ . . . (B 2 ) 


As noted in reference 
the equation 


if 


k= 


11 , this series can be represented by 


AC l 


<y,g 

l+kq 



(B3) 


This result may be interpreted to mean that the successive 
wing-deflection shapes that produce these loads are essentially 
the same, and, hence, this portion of the wing loading may 

1 In reference 11, the constant k was determined by the ratios of the incremental deflec- 
tions 

^ 6] (17) «2 (*?) ^ 

However, for the present method, it will be more convenient to deal directly with the 
incremental lifts. 
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be represented as a single-degree-of-freedom system. Thus 
the incremental coefficient due to this portion of the wing 
deflections may be represented by the expression C L{ q/l-\-kq 
and it will be valid for large positive values of the quantity 
kq even though the series would not converge for positive 
kq greater than 1 . The total lift coefficient for the flexible 
wing may then be written as 

a-^'+TTF^) ,B4 > 

where C Lr is the lift coefficient for the rigid wing. Similar 
equations may be written for the aerodynamic moment 
coefficient. 

Thus it remains to determine C l JC Lr and k through use 
of the influence coefficients. The aerodynamic influence 
coefficients were obtained in the form of the loading coeffi- 
cient G vn at a station v due to a unit angle of attack at 
station n, the angle of attack at the other stations being zero. 
The loading coefficients were obtained by the method of 
reference 10 . By use of this method, aerodynamic influence 
coefficients at four spanwise stations could be found which, 
with proper care, were sufficient to provide desired accuracy 
in determining wing deflections. These coefficients, G vn , were 
obtained by assuming a unit angle of attack at one control 
point, and zero angle of attack at the remaining three, and 
then solving the set of four simultaneous equations consist- 
ing of the a vtl coefficients obtained from reference 10 . For 
instance, G vl can be calculated from the following simulta- 
neous equations 


1 =a u G u +a v2 G2i+a l3 Gsi+a l4 G Al 
0 = a 2 \ G x i -f- a 22 G n + (hz + a 2 t G 41 

0 = cf 31 (7 11 + a32t7 2 i+a3 3 6 r 3 1 +# 34^41 

0 = a 41 Q x ! + a A2 G 2 \ + a 43 f? 3 i + a 44 f ? 41 

where subscripts 1 , 2, 3, and 4 refer to the semispan stations 
77 = 0.924, 0.707, 0.383, and 0 . The resulting lift coefficient 
for each unit angle of attack can then be calculated using 
the equation 

(6 f 4 n+l-848G r 3w + 1.414# 2 n+0.765(7i n ) (B5) 

The experimentally determined structural influence coeffi- 
cients (ref. 12 ) were measured in the form of deflections in 
inches at front and rear spars due to 1 , 000 -pound loads at 
front and rear spars at a number of spanwise stations. 
These were cross-plotted to obtain the influence coefficients 
at the spanwise stations shown in table III. For use with 
the aerodynamic influence coefficients, these coefficients 
were further reduced into the form of a change in streamwise 
angle of attack in radians at station m, due to a 1 , 000 -pound 
load at the quarter-chord position at station r, S mv . 
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The influence deflections due to a unit angle of attack at 
a control station can now be calculated through use of the 
aerodynamic and structural influence coefficients. In order 
to calculate the deflections with sufficient accuracy, an 
integration formula given in reference 10 was used. 

J 1 7T 1 

drj ^ j 2 f i,Vv) sin <f> v 


The total lift for the flexible wing may then be expressed as 


where 


Cl A 


C L =C L \ 1- 


q Cl 


l+kg/ 




k ~ C Li 


(B14) 


where 0„=y— j — ; - and / (17 „) is the value of / ( rj ) at r}=cos rr -- 
l "T l < t1 

For the particular number of control stations used in the 

present case, and since /(t 7 4 ) = 0, and using only the interval 

0<77<1, the integration formula becomes 

f 1 Kv) d v =jz fM L (Be) 

Jo V= 1 

where 

/ 1 =0.1502 
1 2= 0.2 1 76 
73=0.3628 


Thus, the angular deflection at station m due to loading due 
to angle of attack at station n , with q— 1 psi is calculated 
from the equation 


(^0 ) mn 


J> 2 (12) 2 


1000 fri 


3 

2 S m 


G v n I v 


(B7) 


with 


m= 1, 2, 3 
n= 1,2, 3, 4 


The resulting loading at station v due to (e 0 ) mi 
from the equation 

3 

( G \ ) vn ^ ) Gv m ( Cq) mn 
771 = 1 


is calculated 
(B8) 


A similar procedure is used to obtain the aerodynamic 
moment for a flexible wing. 

Wing deflections may also be determined by means of 
wing-deflection influence coefficients in a manner similar to 
that which was used for determining the change in stream- 
wise angle of attack. The deflection at station m due to the 
loading due to angle of attack at station n, with g=l psi is 
calculated from the equation 


with 


fe 2 ( 12) 2 ^ 

J mn 1000 ' 1 v 


rn= 1 , 2, 3 
w=l, 2, 3, 4 


(B 15) 


where Z mv is the deflection in inches at station m due to a 
1,000-pound load at the quarter chord at spanwise station v. 

For an arbitrary angle-of-attack distribution, the deflec- 
tion at station m due to the initial load is 


(^o)m ^ < (^0 )mn&n (Bl6) 

77=1 

The total deflection for the flexible wing may be developed 
in a manner similar to that used for equation (B14). Then 


from which the lift coefficient may be obtained 

7T / I 

(P L jdn = ~ g“ [(^A)4n+l-848(G^) 3n +1.414(G^)2n+0.765((j,i)ln] 

(B9) 

The influence lift coefficient, ( Cl b ) h , due to deflections due 
to the loading coefficient (G A ) yn can be calculated in a 
similar manner. The equation for ( CL B ) n is 

(B10) 

Thus, for a given angle-of-attack distribution, a n , the lift 
resulting from the initial angle of attack and from the first 
and second twist distributions is calculated from the follow- 
ing equations 

4 


( y L R — Xj G Ln Ot n 

77=1 

(Bll) 

4 


r L A =fb (pL A ) n Oi n 

(B12) 

4 



(B13) 




g 

l+kg 


(^ 0 ) m 


(B17) 


in which the k from equation (B14) may be used with suffi- 
cient accuracy. 

APPLICATION TO LATERAL-DIRECTIONAL STABILITY DERIVATIVES 

In much the same manner as in the longitudinal deriva- 
tives, the rolling-moment coefficient resulting from an arbi- 
trary antisymmetric angle-of-attack distribution may be 
expressed as a power series of the dynamic pressure q. 


Gi F —Ci R +Ci A q+Ci B q 2 +Ci c q z + . . . (B18) 

where 


Ci F rolling-moment coefficient for the flexible wing 
Ci R rolling-moment coefficient for the rigid wing 

Ci A q increment in rolling-moment coefficient resulting 

from structural deflections due to the C tR loading 
Ci B q 2 increment in rolling-moment coefficient resulting 

from structural deflections due to the Ci A loading 
For the wing structure considered in the present example, 
the ratios of all terms after the first one, — O lA jCi B , —C 1b /Ci c , 
etc., are essentially equal to a single constant k. The equa- 
tion for rolling-moment coefficient for the flexible wing will 
then be 


‘Ci 




(BIG) 
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The quantities in equation (B19) will now be determined 
through use of aerodynamic and structural influence 
coefficients. 

Antisymmetric aerodynamic influence coefficients were 
obtained through use of reference 18. First, p vn angle-of- 
attack influence coefficients, which represent the angle of 
attack at station v due to a continuous loading function 
having a unit value at station n and zero value at the other 
stations, were obtained from reference 18. These were then 
converted to loading influence coefficients A vn (i. e., the 
loading at station, v , due to a continuous angle-of- 

attack distribution with a unit value at station n and zero 
values at the other stations) by solving for the span loadings 
for a unit angle of attack at one control station. This was 
done for each of the spanwise control stations 1, 2, and 3 
located at 77=0.924, 0.707, and 0.383, respectively. The 
resulting rolling-moment coefficient for a unit angle of 
attack at a control station, n, can then be calculated using 
the following equation, which has the form of equation (15) 
in reference 18. 


[^2n+0.707(^4i n +-4 3n )] (B20) 

n lo 


All aerodynamic loads were assumed to act along the 
quarter-chord line except for loads due to the ailerons. In 
an analogous manner to the symmetrical case, equation 
(B7), the change in streamwise angle of attack at station 
m due to the loading resulting from an angle of attack at 
station n, with q=\ psi, may be expressed as 


b\ 12) 2 * 

\^o)mn 1 ()()() 1 

with 

rn= 1, 2, 3 
n= 1, 2, 3 

where 


(B21) 


A vn aerodynamic influence coefficients for the antisym- 
metric loading function at station v due to an angle 
of attack at station n 

The aerodynamic loading due to ( e 0 ) mn is then obtained from 
the summation 


(A A )pn=^2 A vm {e 0 ) 

m = 1 


mn 


(B22) 
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C, A =S(C^) s a„ 

71=1 

(B25) 

{C, B ) n a n 

71 = 1 

(B26) 


The rolling-moment coefficient for the flexible wing for any 
desired angle-of-attack distribution, a n , is then obtained by 
substituting values from the above equations into equation 
(B19). The reference rolling-moment coefficients with 
the aerodynamic loads at the quarter chord were calculated 
through use of the previous equations (with g=l psi) and 
are tabulated in the following table: 


n 

1 

2 

3 

C, n 

0. 120 

0. 298 

0. 326 

C>, n 

-. 0376 

0098 

-. 0376 


. 00856 

. 01553 

. 00796 


In order to determine the rolling-moment coefficient for 
a specific derivative, the wing angle-of-attack distribution 
must be known. For a unit sideslip angle, the whig con- 
tribution to (\^ may be separated into loadings due to two 
angle-of-attack distributions (ref. 19): (1) a constant span- 
wise angle-of-attack distribution due to the difference in 
effective velocity acting on each wing panel for the yawed 
attitude and (2) an angle-of-attack distribution proportional 
to the wing dihedral angle, which for this case may be 
considered to be due entirely to the upward deflection of 
the wings resulting from the level flight symmetrical air 
loads. 

For the rigid airplane, the contribution to due to the 
differential velocities acting on each wing panel wyis esti- 
mated from reference 20 as C t(} jC L = —0.129. Substituting 

values for a unit angle of attack into equations (B24), (B25), 
(B26), and (Bl9) gives 


C, 

C L 


= -0.129 


1- 


0.1451 

0.744 


, . 0.0321 
1+ 0. 1451^0 


= -0.129 


/ _0 1 195g_\ 

\ l+0.221g/ 


from which the incremental rolling-moment coefficient may 
be obtained. 

(COnha { {A A ) 2n +0.707[(A A ) ln +(A A ) 3n ] } (B23) 

The next incremental rolling-moment coefficient (Ci^) n 
due to deflection resulting from t lie loading coefficients 
(A A ) vn is calculated in a similar manner. 

Thus, the rolling-moment coefficients resulting from the 
initial angle of attack, and from the first and second twist 
distributions can be calculated from the following equations 

n= 1 


The dihedral angle, measured in a plane parallel with the 
YZ plane, was calculated through use of the influence 
coefficients for symmetrical loads and, for the flight range of 
interest, could be expressed as 

^=0.1035g 

The calculated distribution of dihedral angle, normalized 
with respect to the value at the tip, is given in the following 
table: 


V 

0. 383 

0. 707 

0. 924 

r 

. 520 

. 900 

. 992 


(B24) 
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For the wing at a unit angle of sideslip, these values of 
dihedral angle represent changes in the wing angle of attack, 
and the rolling-moment coefficient due to dihedral angle 
may be determined from equations given previously. 


The total value for the wing contribution to becomes 

c L C L + C L 


c 

-^=-(0.1035g) (0.557) 


0.1197 
0.557 q 
i i 0-02664 
1+ 0.1197 


= -0.05772 ( 


0.2142 \ 
1 + 0 . 2222 / 


Since only first-order effects have been considered, the 
changes in dihedral angle due to antisymmetric loads have 
been neglected. 


C 

+=-0.129 

Cl 



0. 195g \ 
l+0.221g/ 


— 0.0577# 


o 


0.214g 
1 +0.222 g 



While the two load distributions yield somewhat different 
values of the ratio CiJCi E , the values of k are essentially the 
same. This was also found to be true for other types of load 
distributions. Note that the values of C ia obtained from 
equation (B27) are not the same as those given in table V 
since only the contribution of the wing has been considered 
here. 


APPENDIX C 


EVALUATION OF LONGITUDINAL STABILITY DERIVATIVES 


Wing lift-curve slope. -This was determined theoretically 
from reference 10 and Appendix B using the section lift-curve 
slope of 5.71 as determined from wind-tunnel data supplied 
by the Boeing Airplane Company. This derivative should 
not be confused with the one in reference 21, which includes 
inertial effects. The present derivative is the one which 
would be evaluated experimentally in a wind tunnel with a 
flexible model. 

Tail lift-curve slope. -This was determined in a manner 
similar to that of the wing, except that fuselage bending was 
included in the flexibility of the tail. It was found that the 
principal reduction in the tail lift-curve slope was caused by 
fuselage bending. Other factors included were the rate of 
change of downwash (ref. 22) and a tail efficiency factor of 
0.95. 

Body and nacelles lift-curve slope. -This was determined 
from wind-tunnel data supplied bv the Boeing Airplane 
Company. 

ft* 

Lag in wing downwash. -Only the tail contribution was 
considered and was determined in a manner similar to that 
described in reference 23, including aeroelastic effects de- 
termined by the method of Appendix B. Although this 
term has small effect from the standpoint of lift, it is im- 
portant in the calculation of C m . 

Normal acceleration. — Since normal acceleration is related 

to a in the equation n= — y ( S—a ), effects of structural 

deflections due to normal acceleration were included in deriva- 
tives in a and 6 by the method of Appendix B. The principal 
contribution to this derivative is from the lift of the wing. 
Contribution of the tail varies according to the amount of 
fuel in the rear main tank and is from 10 to 20 percent of the 
total. 

° L e 

Lift arising from angle-of-attack distribution due to pitch- 
ing velocity (primarily a tail contribution) was determined 
using the method of Appendix B. 
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Normal acceleration. -Same as normal acceleration part 
of C'La} but of opposite sign. 

Lift resulting from angle-of-attack distribution caused by 
structural distortion due to rotational inertial loads was 
determined by the method of Appendix B. This lift is con- 
tributed primarily by the wing but the total effect is small. 



Rigid-airplane value w r as obtained from low-speed wind- 
tunnel data supplied by the Boeing Airplane Company and 
w r as assumed to be constant with Mach number. This is 
justified because the theoretical increase according to the 
Prandtl-Glauert rule is usually compensated by the pressure 
losses at the elevator hinge point (ref. 24). Aeroelastic effects 
of lift and moment on body bending and of lift on stabilizer 
distortion were included. Stabilizer distortion due to elevator 
pitching moment and elevator distortion were neglected. 

Cm. 

The derivative w ; as determined by multiplying the pre- 
ceding C L derivatives by the distance in mean aerodynamic 
chord lengths from their theoretical centers of pressure to the 
center of gravity. The value of C m ^ for body and nacelles 
was obtained by subtracting the theoretical Cm a for wing 
alone from wind-tunnel values of C m for wing, body, and 
nacelles. It should be noted that this also includes the change 
in C vla due to the difference between theoretical and experi- 
mental C„, a of wing alone, a difference which is principally 
due to a somewhat higher loading near the root for the 
experimental than for the theoretical case. Since this addi- 
tional loading occurs near the wing root, it does not affect 
the aeroelastic calculations and, hence, is appropriately 
added to C„ la in the form of a correction for body pitching 
moment. 

(' CCC 

' niQ, ' rrig , m S e 

These moment derivatives were determined by multiplying 
the corresponding lift derivatives by their moment arms in a 
manner similar to that described for C nia . 


APPENDIX D 


PREDICTED LATERAL-DIRECTIONAL AIRPLANE RESPONSE 

EQUATIONS OF MOTION 

The three lateral-directional equations of motion with 
respect to stability axes given in Appendix C of reference 3 
for a rigid airplane will be modified foi use for a flexible 
airplane. For the flexible case, additional terms must be 
added to take account of aerodynamic forces resulting from 
structural deflections due to inertial loads C v ^ C n C Y ^ 
and C Y 0 and dead-weight loads (C Yif> ). A more complete 
evaluation of these additional derivatives would also have 
included other derivatives dependent on Ip, \ }, \p-\-p, and <p. 
However, by the substitution of typical numerical values in 
the transfer -function coefficients, the cumulative effects of 
these other derivatives on the transfer-function coefficients 
were found to be negligible even though some of them were 
of appreciable size (e. g., (C n ^) and C n - ) were several times 
larger than C nr {b/2V ) for the largest values of M). The 
rolling moment, yawing moment, and side-force equations, 
with flight-path angle assumed zero, may be expressed as 

(l x D*-q 0 SbC, 9 D 2 -q 0 SbC lp 2 b v D) ^+(-/v Z />- 

g„SbCi.D 2 —q 0 SbCi r A D\-q 0 SbC^=q 0 SbC ls d (Dl) 

(~I xz D 2 ~q 0 SbC np A Dy+(l y D 2 -q a SbC ni D 2 - 

q 0 Sbr„ r A /;) i-q 0 SbC n Ji= io SbC n l> (D2) 

(-W-q 0 SC Yv )<p+ ( mVD-q 0 SC r .D ) 

(mVD- q 0 SC y .D- q 0 SC Yfi )fi= q 0 SCy s 8 (D3) 

The remaining inertial deflection derivatives were then 
combined with the mass terms to form effective inertias for 
the flexible airplane. Thus 

Ix F — lx (J 0 Sb C\ „ 

h F =Iz q 0 SbC n 

Qo ^ sy 

C Y . 

V 0 

q 0 S n 

or m—~ C Y - 
V i 

Since the derivatives C Y - and C Y • each represents aerodynamic 
forces due to lateral acceleration, they are of equal magnitude; 
and hence m F may be expressed in terms of either derivative- 
The three equations (Dl), (D2), and (D3) can then be 
written in a more convenient form by dividing them by 
and m F V, respectively, and introducing new 

symbols 


(D 2 -L p D)cp+(-r XF D 2 -L r D)^-L^=Ls8 (D4) 

{-r ZF n 2 -N p D)<p+ (D 2 -N r B)t-Ntf=N s 8 (D5) 

— 1 &8-{- (D — 1 &)P = ^ (D6) 

Note that the term Y<p (which equals g/V) remains the same 
as for the rigid case since the derivative C Yip is due to dead- 
weight loads that are distributed in the same manner as the 
lateral acceleration loads. 

TRANSFER FUNCTIONS 

Transfer-function coefficients for <p/8, \///8, and fi/8 can be 
obtained by simultaneous solution of the three equations 
(D4), (D5), and (D6). 

Denominator coefficients. The denominator of the trans- 
fer functions, A, may be expressed as 

A = D(C a D*+ CJP + CJP+CiD+Co) 

where 

— 1 V X F r Z F 

C z =-L p '-N r '-Y fi (l-r XF r ZF ) 

C 2 =Ne'+ Y p (L/+N/)+L p N-N p L r 
Ci = LpNp — NpLp— Y v W+Y p (L r N p -N r L p ) 

C 0 = ) <p(LpN r — N&L r ) 

For the moderate angle-of-attack range considered, r Xp and 
r Zp are small quantities and hence C±~ 1. The denominator 
can then be expressed in factored form as 

A=D(D + D s )(D+D r )(D 2 +c 1 D+c 2 ) 

where D s and I) T represent the spiral and rolling modes, 
respectively, and c x and c 2 are coefficients that define the 
damping and frequency of the Dutch roll mode. 

Numerator coefficients. — With p=D<p , and r=D\p, the 
numerator coefficients of the following transfer functions can 
be expressed in terms of the derivatives. 

p_ D(AJP+A 2 D 2 +AiD) 

8 A 

where 

A 3 =L/ 

A 2 = — ) pLs' -j-JVsLr — LsN r -j- ) &Lp 
Ai=L&N0 AVZ'/sT 5 & (LdN r — A oLr) F oC^rA $ — A T L$) 
r_D(B 3 D*+B 2 D 2 +BiD+B 0 ) 

8 A 
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where 

B 3 = AY 

B 2 =L^N v — NiL v — ) oAY + } sNb' 

B x = Y0(NiL 9 —LiN p ) + Yi(LfiNp-NfiL p ) 

7 ? 0 = Y^LsNfi—NsLp) 

P D(E z D z +E 2 D 2 +E l D+E 0 
8 A 

where 

Ez = Y& (1 — TxfTzf) 

E 2 =-Y s (N r '+L p ')-N 5 ' 

Ei = Y JLi + N&L P — Lb A p + I s(L p N r — A v L r ) 

£'o= Y^NsLr—LiNr) 

Simplifications of the transfer functions. — When transfer 
functions are evaluated from measured data, it is desirable 
to use as simple a form as possible which will still adequately 
fit, the data. Calculations indicated that the spiral mode 
factor D s was very small and could be neglected for the 
frequency range of interest. Thus the rolling response to 
aileron was simplified to 

p_ A3 (D 2 -\-a\D-\-a 2 ) 

8a~(D+D T ) (Dt + CyD+Ct) 

where 


A 1 

a2 ~A 3 

To determine estimated values for use with the curve fitting 
of the measured responses on the analog computer, a partial 
fraction form of the transfer function was advantageous 
to use. 

j) Gp H p D-[-Jp 

d^D+D^D'+dD+Ct 

In obtaining approximate relations for r5 r and P/8 r , it was 
desirable first to write the transfer functions in partial 
fraction form. 

r F r G r | H r D-\-J r 

8 r D-\-D s D+D r D 2 +CiD-\-c 2 

P Ffi G$ . HpDYJp 

8 t D-\-D s D-\-D r D 2J rC\D-\-C2 

By the substitution of typical numerical values, all numerator 
terms were found to be negligible except H r and J^. Thus, 
the transfer-function coefficients could be simplified to 

r = H T D 

8 r D 2 -\-CiD-\-C 2 


Transfer functions for the rear part of the fuselage and 
the wing tip. — In addition to the transfer-function coeffi- 
cients for quantities at the center of gravity, the coefficients 
for r t /8, p t /8, and n wt /8 are needed. The evaluation of these 
quantities includes responses at the center of gravity and 
also responses due to pseudostatic structural deflections at 
the particular location. The equation for yaw rate at the 
rear part of the fuselage is 

r t _r . DA\pt 

V 

where A\p t represents the total change in angle of yaw in 
radians at the rear part of the fuselage due to structural 
deformations resulting from a control input. The following 
quantities were included in determining A \l/ t /8 

A^_A^ P . Ai 1/ r Ajt $ . A^ f 
8 p <5 ' r 5" 1 " y 8^ 8 

For instance A \ptlp represents the pseudostatic change in 
angle of yaw at the rear part of the fuselage due to structural 
deformations resulting from a unit change in p. 

Similarly, the equation for roll rate at the rear part of 
the fuselage is 

Pt_JP \ Dk<Pt 

8 8 ^ 8 

where A <p t represents the total change in angle of roll at 
the rear part of the fuselage due to structural deformations 
resulting from a control input and may be expressed as 

A (p t A*, P . r . A <p t A iPt 

8 “ P 8^~ r 8^ f 5 5 

The equation for acceleration at the left wing tip is 
w* wi __ L)p . D 2 z wt 

~T~ 126” 

Calculations indicated that the effect of D 2 z wt /8 was small 
for the frequency range up to the first antisymmetric bending 
frequency and, hence, was neglected. 

Conversion of transfer-function coefficients from stability 
axes to body axes. — In order to compare predicted results 
with measured frequency responses which were obtained 
with respect to body axes, the predicted transfer functions 
were converted from stability-axis to body-axis form. Since, 
for the angle-of-attack range considered, the quantity 
was essentially unity, only the numerator terms of the 
transfer function need be modified. The equations for con- 
version from stability axes to body axes are (e. g., ref. 3). 

p B = p cos a — r sin a 

r B =r cos a-\-p sin # 


P = Js 
8 D 2 -\- C\E-\-c 2 


tan p B = 


tan p 
cos a 
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For small angles of attack, the equations may be simplified 
to 

p B =p—ra 

r B =r+pa 

These relations will also hold for derivatives of the angles. 
Hence, the numerator coefficients for p/8 and r/8 may be 
converted as follows: 

A iB =Ai— B T a 
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